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Abstract 

We consider the algebraic Bethe ansatz solution of the integrable and isotropic XXX- 
5 Heisenberg chain with non-diagonal open boundaries. We show that the corresponding 
-RT-matrices are similar to diagonal matrices with the help of suitable transformations 
independent of the spectral parameter. When the boundary parameters satisfy certain 
constraints we are able to formulate the diagonalization of the associated double-row 
transfer matrix by means of the quantum inverse scattering method. This allows us 
to derive explicit expressions for the eigenvalues and the corresponding Bethe ansatz 
equations. We also present evidences that the eigenvectors can be build up in terms of 
multiparticle states for arbitrary S. 
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1 Introduction 



The possibility of constructing SU(2) invariant Heisenberg chain with arbitrary spin-S* solvable 
by Bethe ansatz methods was a remarkable achievement of the representation theory underlying 
the associative algebra describing the dynamical symmetry of quantum integrable systems pQ. 
It turns out that the Hamiltonian of such spin-5 1 XXX Heisenberg magnet [2j commutes with 
the transfer matrix Tg(A) of a 25* + 1 state vertex on the square L x L lattice PJ Ej- This 
connection is based on well known relationships between one-dimensional quantum spin chains 
and two-dimensional statistical mechanics models whose Boltzmann weights satisfy the Yang- 
Baxter equation [3J Ej . 

The row-to-row transfer matrix Ts(A) of such 2S+1 state vertex model can be conveniently 
written as the trace, over an auxiliary space A = C +1 , of an ordered product of Boltzmann 
weights. More specifically, 

Tg{\) = Tr A [Tf(\)\ Tf (A) = £S(A)£2_ 1 (A) . . .£<g(A), (1) 

where A is the spectral parameter and A represents the horizontal degrees of freedom of the 
vertex model. 

The Boltzmann weight C ab (A) is solution of the Yang-Baxter equation 

4?(A - /i)Tj s )(A)7f V) = ^(aO^C*)^ - /i), (2) 

invariant relative to the SU(2) Lie algebra. It can be viewed as (2S + 1) x (2S + 1) matrix on 
the auxiliary space whose elements are operators acting non-trivially only on the 6-th factor of 

L 

the Hilbert space Its explicit expression in terms of the spin-S* SU(2) generators 

6=1 

S a = (ss,sisz) is m U El 

££ » w=( A +2 , 5) ^n^|n^^#_-, (3) 

1=0 k=l + l 1 n=0 ' n 

where xi = + 1) — S(S + 1) and r\ is the so-called quasi-classical parameter. 
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Besides the Yang-Baxter equation the operator C^iX) satisfies other relevant properties 
such as 

Unitarity: C$(\)£$(-\) = Cs(A)Id ® Id; (4) 
Parity invariance: P^C^' \X)Pi 2 = 4.2p(A); (5) 
Temporal invariance: 4f ) (A)' 1 ' 2 = 4? (A); (6) 

Crossing symmetry: 4?(A) = (-1) 25 N V C[f(-X - r/)' 2 (7) 

- ss{-X-ri) 

25-1 

where functions Cs(A) = (2Sr]) 2 -X 2 and q s (X) = ] \ (X+kr]). Here Id is the (25+1) x (25+1) 

k=l 

identity matrix, Pyi is the permutation operator, t a denotes transposition on the a-th space, 

1 2 

V— V ® Id and V— Id <8> V. The matrix V is anti-diagonal whose non-null elements are 
V M = -(-l)%2S+2-r 

This notion of integrability has been extended to include integrable open boundary condi- 
tions [HllZj- In addition to the Yang-Baxter solution C ab (A) determining the dynamics of the 
bulk one has to introduce (2S + 1) x (25* + 1) f^-matrices K$(X) whose elements represent 
the interactions at the left and right ends of the open spin chain. Compatibility with bulk 
integrability demands that these matrices should satisfy the reflection equation given by [7j 

4? (A - /i) Ks (A)4?(A + fi) Ks (ji) =K S {^{X + //) Ks (A)4?(A - //), (8) 

where Ks (A) = K S (X) ® Id and Ks (A) = Id ® K S (X). 

In the case of open boundaries the analogue of the transfer matrix is the following double- 
row operator 



ts{X) = Ti A 



K^(X)T^(X)K { i\x) TH-A) 



+5), 



(9) 



where Kg \X) can be chosen as one of the solutions of the reflection equation (jHJ). The 
other matrix Kg ■ (A) can be directly obtained from K s (A) thanks to the extra relations ^ 
EJ) satisfied by the operator C ab (A). Following a scheme devised in ref.jS] this isomorphism 
becomes 

K s + \X)=\K S '\-X- V )] t . (10) 



The understanding of the physical properties of the XXX-S" open chain includes necessarily 
the exact diagonalization of the double-row operator (|5)l. If the i^-matrices are diagonal 
this problem can be tackled, for example, by an extension of the quantum inverse scattering 
method [7j and the use of fusion hierarchy procedures EI]. The same does not occur 
when the i^-matrices are non-diagonal due to an apparent lack of simple reference states to 
start Bethe ansatz analysis. In spite of this difficulty, progresses have recently been made 
for the anisotropic version of the S = \ Heisenberg magnet usually denominated the XXZ 
spin chain. These achievements have been made either by a functional Bethe ansatz analysis 
[TT] or by means of the algebraic Bethe ansatz method The latter approach has been 

based on earlier ideas developed in the context of the eight vertex model ^3] • in particular, it 
was argued that the spectrum of the open XXZ chain can be parameterized by Bethe ansatz 
equation provided certain constraint between the parameters of the Hamiltonian is satisfied. 
Part of the conclusions were achieved with the help of a numerical study of the spectrum for 
finite values of L ^3]. More recently, new results have been obtained in ref.[15j by exploring the 
description of the open XXZ spin chain in terms of the Temperley-Lieb algebra. The extension 
of all such analysis for integrable Heisenberg chains with arbitrary spin-S" appears to be highly 
non-trivial and it is indeed an interesting open problem in the field of integrable models. 

In this paper we would like to take some steps towards the direction of solving the isotropic 
higher spin Heisenberg model (jUJ) with non-diagonal open boundaries. We show that the 
double-row transfer matrix operator associated to the integrable XXX-S* Heisenberg chain can 
be diagonalized by Bethe ansatz at least when the respective i^-matrices parameters satisfy one 
out of two possible types of constraints. We find that the roots of the Bethe ansatz equations 
are fixed by integers n < 2SL that play the role of standard particle number sectors. This 
feature shows that the Hilbert space has a multiparticle structure which should be useful to 
determine the nature of the ground state and excitations unambiguously. 

The outline of this paper is as follows. In section |21 we argue that the non-diagonal K- 
matrices of the XXX-5* Heisenberg model are diagonalizable by spectral independent similarity 
transformations. In section |3] suitable quantum space transformations are used to show that 
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the diagonalization of is (A) is similar to an eigenvalue problem with diagonal and triangular 
K- matrices provided that certain constraints are satisfied. In section 0] we discuss the quantum 
inverse scattering method for the latter system, presenting the corresponding eigenvalues and 
Bethe ansatz equations. Explicit expressions for the eigenvectors in terms of similarity trans- 
formation acting on creation fields can be written for spin | and 1. In section El our conclusions 
and further perspectives are discussed. In Appendix A we summarize certain properties of the 
X-matrices. In Appendices B and C we discuss the one and two particle analysis of the eigen- 
spectrum as well as auxiliary expressions for S = |, respectively. Finally, in Appendix D we 
exhibit general relations concerning the one-particle unwanted terms and the two-particle state 
construction for arbitrary S. 



2 The i^-matrices properties 

The most general reflection if-matrix associated to the open XXX-S" Heisenberg chain possesses 
three free parameters. For S — | it is given by ^B] 



(A) 

2 



while the isomorphism (jlUj) implies that 



K[ + \\) 

2 



-cL 



V 



1 — - 
* + l 



(ii) 



e + + l+t 



(12) 



where £±, c± and d± are six free parameters. 

A remarkable characteristic of these iiT-matrices is that they can be diagonalized by sim- 
ilarity transformations which are independent of the spectral parameter A. More precisely, it 
is possible to rewrite the equations (jlipi2j) as 



K^\\) = P \yg\ 



(-)r.(-) 







(-) 



-i 



(13) 
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and 



o e+ + 1 + * 



(+) 



(14) 



where Qg refer to appropriate (2S + 1) x (2S + 1) matrices. In what follows we will represent 
them in terms of the standard Weyl basis by the expression 

2S+1 



Q {±) 
y s 



(±) 6 



(15) 



In the specific case of S = h the expressions relating the off-diagonal and the diagonal 



elements of g[ are 



1 + e ±V 1 + c±d± 



1 + e ±v 1 + c ± d ± 



(16) 



where e + 
given by 



e_ = ±1. The other variables £± and pi entering in the formulae (|lH|14jl are 

2 

6 ±^± . .(±) 



± 



-e ± yT+ c±d ± . 



(17) 



a/1 + c±d± 

The i^-matrices for S > ~ can be computed either by brute force analysis of the reflection 
equation J7j or constructed by the so-called fusion procedure [18J. Their matrix elements 
expressions become very cumbersome as one increases the value of the spin and this fact has 
been exemplified in Appendix A for spin 1 and | cases. It turns out, however, that we have 
found out that such i^-matrices can be rewritten in a rather compact and illuminating form 
with the help of appropriate spectral independent similarity transformations, namely 



(±)r>(±)n(±)/ 



y s 



18) 



where the overall normalizing constant is p s 



(±) 



,2S 



a/1 + c±d± 



25 



2S 



The diagonal matrix D^\X) is defined by 

25+1 



(19) 



where the corresponding diagonal entries are 

25 

fL~\S; A,£_) = n 
0=1 

2S 

n 

0=1 



L + S + ^-P- sign (a-^-P 



/i +) (S;A,£+) 



C+ + S+ l --P + sign 1-/3^(^ + 1) 



(20a) 



(20b) 



Interesting enough, the novel parameters £± encode both the dependence on the spin value 
and on the variables describing the off-diagonal i^-matrices elements. Specifically, we have 
found that 

4 = ~^=. (21) 



a/1 + c±d± 

Finally, the four elements , , ^ l an d 022^ °f @s are related by the expressions 



9® 



l) 2S V2S(l + e±^/l + c ± d± 



:-i) 2s V2 

VSc± 



(S - l)c±d± + (2S -!)(! + e ±v /l + c ± rf±) 



1 + e ± yl + c±ci ± 

and the remaining elements are obtained by the following recurrence relations 

(±) y/2S(m - 1)(2S + 2 - ggg^g^ - y/2S(l - 1)(2S + 2 - jjgjgggL 
^ 2S(m-i)«$ 
for I ^ m = 1, . . . ,2S + 1 while for m = / we have 

ffSff&iw ^2(25 - 1)(Z - 2)(2S + 3 - t)g&(S)g}% + 2(1 - 2)g^g\% 



a {±) 
yi,i 



a {±) 
yi,i 



(22) 
(23) 

(24) 

(25) 



v /25(25 + 2-/)(/-l) y i ± 1 ) 
An important feature of our construction is that the matrices are itself representations, 
without spectral parameter, of the monodromy matrix associated to the Yang-Baxter algebra 

( Q\ —————————— 

(J2J) generated by the operators C ab (A). In fact, the matrix (|15|24|25jl with four free parameters 
are the widest possible class of non-diagonal twisted boundary conditions compatible with 
integrability for the XXX-S* spin chain |T5] . An immediate consequence of this symmetry is 
the commutation relation 



4?(A),S 



which will be of great use in next section. 



(±) 

5 



0. 



(26) 



3 The eigenvalue problem 



The purpose of this section is to show that the eigenvalue problem for the double-row transfer 
matrix operator ts{X), 

t g {\) |V>) = As(A) \if>) , (27) 

associated to the XXX-S* chain with two general non-diagonal open boundaries can be trans- 
formed into a similar problem with only one genuine non-diagonal i^-matrix. 

In order to demonstrate that we use the decomposition property for the K^~\X) matrix 
described in section |21 and the operator ts(A) becomes 



ts(X) 
Ps 



Tr.4 



(A) 



a (+y 



(S), 



A 



-A) 



(2* 



in between all the fun- 



We now proceed by inserting identity terms of type 
damental operators that appear in the trace (J28|) . By using the invariance of the trace under 
cyclic permutation one can rewrite Eq. (j28j) as 



«s(A) 



Tr 



A 



d£\\)tx»(\)kP(\) 



where Tf } (A) = £^(A)£^/_ X (A) . . . C^((X). The new operator £% j (A) and fT-matrix Kf'(X) 
are given in terms of unitary transformations acting on the auxiliary space by the expressions, 



Tf(-A) 



(29) 



; (-), 



and 



!(+) 



-1 



(+) 



4 _) (a)<? 



(+)_ 



(30) 



(31) 



It turns out, however, that the gauge transformation (|3U|) on the £^-(A) operators can be 
reversed with the help of a second transformation on the quantum space [19 . In fact, one can 
use property (j26|l to define quantum space matrices Vj acting non-trivially only at the j-th site 



Id 



Id 



rth 



Did 



Id, 



(32) 



such that they are able to undo the transformation (|30|) . namely 

V^5J](X)Vj^C^(X). (33) 
We now can use this property in order to define a new double-row transfer matrix operator 

L L 

having only one non-diagonal if-matrix 

(35) 

Clearly, the operators tsr(A) and is (A) have the same eigenvalues while their eigenstates |^) 
and are related by a similarity transformation, 

L 

\^)=U V A^)- (36) 
i=i 

Though this framework clearly brings a considerable simplification in the original eigenvalue 
problem, it is not enough to make the diagonalization of the double-row operator is (A) with six 
free boundary parameters amenable to a standard Bethe ansatz analysis. This is because the 
i^-matrix K s is generally non-diagonal which still imposes us the difficulty of finding suitable 
reference states needed to begin the Bethe ansatz computations. However, a great advantage 
of this formulation is that one can easily identify the existence of at least three cases of physical 
interest in which the standard SU(2) highest weight states could be used as pseudovacuums 
to build up the whole Hilbert space. The simplest occurs when one of the boundaries is free, 
say K S ~'(X) = Id while the other is still arbitrary with three free parameters. The next 
one is when the .fT-matrices K s ± \\) are diagonalizable in the same basis, i.e. Q s + ^ = Q s ~^ 
which implies that we have altogether four distinct couplings say c_, oL and £±. This includes 
the important symmetric situation where the left and right i^-matrices are the same but 
arbitrarily non-diagonal. As far as the Bethe ansatz technicalities are concerned the most 
general case in which SU(2) highest weight vectors can be used as a reference state is when 
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ts(X) 



Tr 



A 



D 



{ s + \X)Ti s \X)K ( s -\x) 



(S), 



A 



-A) 



the effective K s \X) i^-matrix becomes either upper or lower triangular. This leads us to an 
open integrable system with five free couplings since such condition imposes certain constraint 
between the parameters c± and d±. Substituting the representation (JTSj) in Eq. (}3*T]) and after 
some algebra we find that there are two possible classes of restrictions satisfying the above 
mentioned triangularity property It turns out that these constraints depend only upon the 
variables c±, d± and their expressions are, 

1 + 6,^1 + c_gL = 1 + e +v /l + c+d+ 
c_ c + 



(ii) ~ d - l+e+Vl+cA ^ 

l + e_ A /l + c_rf_ c+ 



Depending on the ratio e = — the zeros entries of K S ~\X) are either bellow or above the 
principal diagonal. This feature has been summarized in Tabled for each manifold. Note that 
the diagonal elements of the triangular matrix Kg (A) will necessarily be the eigenvalues of 
K s (X). By considering decomposition (fTH|) we conclude that such eigenvalues are exactly the 
entries of the diagonal matrix pgDg(X). 

Considering the above discussions, we find that the formulation of a Bethe ansatz solution 
for the eigenspectrum of is(X) on the parameters manifolds (I) and (II) is certainly worthwhile 
to pursue. It will leads us to benefit from the knowledge of the exact spectrum with five out 
of six possible boundary couplings, a considerable number of free parameters at our disposal. 
A fundamental ingredient in the algebraic Bethe ansatz is the quadratic relations satisfied by 
the matrix elements of the double-row monodromy matrix defined by [7j 



TfW = Tf(X)Kir\\) [TX"(-X)\ . (39) 
and consequently the double- row operator ?s(A) can be written in the form 



-(5) 



-1 



„(+) " LlA 
Ps 



Df\X)Tf{X) . (40) 



Taking into account the property f!26l) we see that the effective Kg (A) matrix satisfies the 
same reflection equation (jHJ as the original i^-matrix Kg (A). As a consequence of that and 



the fact the entries of Kg~\\) are c-numbers it follows that X^f (A) is also a solution of the 
reflection equation, namely 

1 2 2 1 

£g>(A - a*) (A)4?(A + a*) (/x) =7<? (Ai)4?(A + a*) (A)4?(A - A*)- (41) 



In the next section we will explore such quadratic algebra together with the existence of a 
pseudovacuum on which T A (A) acts triangularly to present the expressions for eigenvalues of 
is (A) as well as the corresponding Bethe ansatz equations. 



4 Algebraic Bethe ansatz 

In the next subsections we will consider the diagonalization of the operator is (A) in the most 
general restrictive condition (I) or (II) by an algebraic formulation of the Bethe ansatz. The 
other two situations mentioned in section El are special cases and the corresponding eigenvalues 
and Bethe ansatz results can be derived from the results, for example, obtained for manifold 
(I). This is obvious when Gs^ = Gs^ and in the case K s \X) = Id one needs to consider 
the limit — > oo with fixed p s — 1 in the results to be given bellow. 



4.1 The spin-7> solution 

Here we shall consider the diagonalization of the double-row transfer matrix ii(A) by means 
of the quantum inverse scattering method 0IZI- The corresponding bulk Boltzmann weights 
(J3J) are those of the isotropic six-vertex model, 



4l } (A) 



X + 7] 
A 77 

T] X 









(42) 



y X + r] J 
Following the remarks of section El we are assuming that the boundary couplings c± and 



d± satisfy one of the two possible constraints described by Eqs.flc 



In this situation 
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the effective Kg \X) 7^-matrix is triangular and its diagonal entries are proportional to the 



eigenvalues /• A, £_). Without losing generality one can clearly consider the case in which 
k[ \X) is upper triangular, and after some simplifications in Eq. (|31|) we find that 

2 

K ( f\X)=ep ( [ 



o 



(43) 



The off-diagonal term in Eq. (j43j) is not expected to affect the eigenvalues of ii(X) but it 
will certainly be relevant in the structure of the eigenvectors. The explicit expression for a 12 
has been presented in Appendix A. As discussed in section Ola direct consequence of the upper 
triangular property of Ks{\) is that the following SU{2) highest state vector 

L 



|0s> 



i=i 



\S,S) 







V / 



(44) 



25+1 



is an exact eigenvector of the double transfer matrix ts(X). 



Oi ) can be used as pseudovacuum to build up the other eigen- 



This means that the state 

vectors of ii(X) following the strategy of the algebraic Bethe ansatz approach |3J[7]. A main 

2 — ( l ) 

step in this method involves writing the double-row monodromy matrix T J (A) in the 2x2 

form 

' A(\) B(X) 
C(A) D(X) 

By using the intertwining relation (j41j) and following the procedure devised first by Sklyanin 
[Zj one can derive the commutation rules 



7p(V2) 

1 A 



(45) 



A{X)B(n) 



D{X)B(ji) 



Qu-A + ry)(/x + A) 



2ll 



+ 



[fi + X + r]) 

(X + Li + 2rj)(X- li + t]) 

(A - /i)(A + fl + T}) 

4?7/i(A + 77) 
(X + Li + r])(2X + t])(2li + i]) 
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B(fi)D(X) 



(2/x + Ty) (/i-A) 
2ri(X + r]) 



(46) 



(2X + V )(X-li) 



B{X)D{p) 



B{\)A{n), 



(47) 



where the new field -D(A) is introduced in order to simplify the commutation relations. It is 
given by the following combination between the operators A(X) and -D(A), 

V 



D(X) = D(X) 



-A(X). 



(4£ 



2X + 7] 

In terms of the operators A(X) and D(X) the double-row transfer matrix eigenvalue problem 
can now be written as 



A(X)\^) + ft\^X^ + )D(X)\^) 



1 



P 



(+) 



(49) 



while the action of the fields A{X), D(X) and C(X) on the reference state Oiy are given by 

L 





A(X) 

5(A) 
(7(A) 



2 2 



(A + ?y) 2 



0. 



(-) 



Ci(A) 



/i H (^;A,^- 



2X + r] 



A 2 



a (a) 



The fields -B(A) are interpreted as a kind of creation operators over the pseudovacuum 
and the multiparticle Bethe states |Vv(Ai, . . . , A n )) are supposed to be given by 



(50a) 

(50b) 

(50c) 




\$ n (X u ...,X n )) = B(X 1 )...B(X n ) 



0i 



(51) 



The rapidities Xj will be determined by solving the eigenvalue problem with the above 
ansatz for the eigenvectors. This is done with the help of the commutation relations (J4bl 
W7\ to move A(X) and D(X) in Eq. (J49|) over the creation fields until they reach the reference 
state Oiy. The terms proportional to the eigenvectors (|51|) are easily collected by keeping 
only the first pieces of the commutation rules. After using expressions (j50aH50b|) and some 
simplifications we find that the final result for the eigenvalues are 



Ai(A) 

2 2 



+ 



(A + yf 

a (A) 

A 2 



a (A) 



£ 2(A + 7y)(A + £7ye-)(-A + ^+ 
(2A + r])i] 2 

e2A(-A + erj^ - t?)(A + rjC+ + rj) 
(2A + r])r] 2 



" (A i -A + 2)(A i + A- 



n 



L(A,-A-2) (A, + A + l) 



+ 



fV ^.-A-f) (Aj + A + f) 
fi (Ai-A-I) (A i + A + D l J 
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where we have used the values of /• (§; A,£±) taken from Eqs. (j2()a[ I20bj) and performed the 



V 



Xi — — on the rapidities. 



A n )) can be canceled out 



displacements A; 

The remaining terms that are not proportional to IV^Ai, 
by imposing further restrictions on the rapidities Xj. These are known as the Bethe ansatz 
equations which in our case are given by 





~ 2 


2L ^ 




V 

2 - 





Xj+ £7?e--i 



-q (A, - A^ + 77) (Aj + Aj + 77) ( . . , , 



Aj - vt+ ~ 1 / if ( A i — Ai — 77) (Aj + A* - 77) 



We now can derive similar results for the open spin-1 chain that commutes with the double- 
row transfer matrix ti(A). The corresponding Hamiltonian is proportional to the first-order 
expansion of 1 1 (A) in the spectral parameter pj HHj 

L-l 



V 



i=i 



(54) 

where 0^, 0^ and cr* are the Pauli matrices with = \ (a* ±ia^). Its eigenvalues E n (X) 
are obtained in terms of the rapidities Xj that satisfy the Bethe ansatz equation (|53J1 by the 
following expression 



E n (X)=2 V J2 



X 2 _ 221 

fe=l A fc 4 



L 1 

77 77 



P 

1+ 1 



(+) 



(-) 



(55) 



We would like to close this section with the following remark. The ferromagnetic 77 < 
Hamiltonian (J53j) is known to describe the stochastic dynamics of symmetric hopping of parti- 
cles in one dimension provided that certain relations are satisfied by the boundary parameters 
PU] . More specifically, letting (2(7) be the rate of injection (ejection) of particles at the left 
boundary and 5 (ft) the corresponding rate at the right boundary we have [2lH |2*T| 



and 



a — 7 



(3-5 



Oi 



d- 
2^ 



c+ 
'277^ 



7 



d + 



(56) 



(57) 
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The above particular parameterization of the boundary parameters c± and d± satisfies the 
constraints (I) or (II) for arbitrary values of the particle injection and ejection rates. Though 
the spectrum at this special case have been determined before |2()| l2"T] not much is known 
about the behavior of the wave functions. This information can now be in principle extracted 
by combining the unitary transformation (jH2l ESI) with the multiparticle state structure (|51|). 
This knowledge of the eigenvectors can be used to calculate correlation functions, thanks to 
recent developments made in the quantum inverse scattering method [22 EH] which allows 
us to reconstruct local spin operators in terms of the monodromy matrix fields. We hope to 
return to this problem since this could provide us with new insights on the physics of stochastic 
dynamics of interacting particle systems. 



4.2 The spin-1 solution 



The statistical system associated to the integrable XXX-Heisenberg model with spin-1 is a 
three-state vertex model with nineteen non-null Boltzmann weights given by 

\ 



4S(A) 



o(A) 





























5(A) 





h(X) 























e(A) 





d(X) 





5(A) 











h(X) 





6(A) 























d(X) 





<?(A) 





d{X) 























5(A) 





h(X) 











c(A) 





d(X) 





e(A) 























h{X) 





5(A) 





























a(A) 



with 



a(A) 



A + 277, 6(A) = A, 
X(X- V ) 



c(A) 



2 V 2 



d(X) 



2 V X 



e(A) 



X + r] 



X + T] 1 X + T]' 

g(X)=b(X)+c(X), h{X) = 2r 1 . 



(5? 



(59) 
(60) 
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As before we can consider the situation when the effective i^i(A) matrix is upper triangular. 
In this case, carrying on few algebraic simplifications in Eq.(|3ip we find that 

/ 2 (l;A,^_) 
o 



MX) 



\ 



/i(i;A,^- 
o 





«12(| - ^ 



^23 



2^13 (C 
o <+) * 



A 1 
r, I 2 



«12(| - f ) + K 23 



/a(l;A,ee-) 



(61) 



where the off-diagonal coefficients K12, K13 and K23 have been collected in Appendix A. 

At this point we need to start introducing suitable notation for the double monodromy oper- 
ator T A (A). Here we shall use a representation which can be easily extended to accommodate 
arbitrary spin-S 1 case, 

( A 1 (X) B 12 (X) B 13 (X) N 
?aW= C 21 (A) A 2 (X) B 23 (X) ■ (62) 
v On(A) C 32 (A) A 3 {X) j 

The next step is to rewrite the eigenvalue problem in terms of the double monodromy 
matrix elements. To perform this task is convenient to define new diagonal operators Ai(X) in 
terms of appropriate linear combinations of the fields Ai(X) 0E1!- This is done in such way 
that the action of the new fields on the state |0i) will be proportional to a single bulk term. 
Keeping in mind possible extension to general values of the spin we define, 



A,(X) = A 1 (X), 

A 2 (X) = A^ + ^MX), 

*< A > = l3 ^ + S^ Il(A) + l§AT l2(A) ' 

where the functions /ii(A) and h 2 (X) are the following determinants 



MA) 



a(X) c(A) 
h(X) h{X) 



a(A) h(X) 
h(X) g(X) 



(63) 
(64) 

(65) 



(66) 
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Taking into account the representation (|62|) and the above redefinitions of the diagonal 
fields, the diagonalization of the doubled transfer matrix ti(A) becomes equivalent to the 
problem 



^^ (+) (A)^(A)|^ n (A 1 ,...,A n )> 



Ai(A) , 7 n 



i=l 



A„)> 



with 



^S +) (A) = /r(i;A,e+) + 



(+)/ 



fc(2A) .(+) 



a(2A) 



/r ; (i;A,e+) + 



c(2A) 
a(2A) 



/ 3 (+) (l;A,en 



4 +) (A) = / 2 w (l;A,e + ) + 



(+)/ 



/n(2A) 
R2A) 



/i +) (i;A,e + ), 



^ +) (A) = /ni;A,£ + ) 



(+)/ 



(67) 

(68) 

(69) 
(70) 



Another important ingredient is to determine the action of the double monodromy matrix 
elements on the pseudovacuum |0i). This can be done with the help of the Yang-Baxter 
algebra |2H an d the triangularity properties of both £^(A) and K[~\X) operators upon 
|0i). Following ref.|21| and taking into account Eq. (J61|) we have 



with 



co{-\X) 
^ } (A) 

Jf\X) 



M*) |0i> 
MX) |5i> 
MX) \o t ) 

C 2 i(A) |0i> 

A H (i;A,£f- 

/ a (_) (l;A,e£. 



P^M^A) 

^2~ } (A) 
p^u^A) 



q(A) 2 

LG(A) 

Ci(A) 
e(A) 2 



|0x> 



|0i) 



|0i) 



LCi(A). 

C 3 i(A) |0i> = (732(A) |0i> = 0, 



M2A) , ( _) 



a(2A) 



/r(l;A,e£- 



^i(2A) , H 



where the new function /13(A) 



/i 2 (2A) 

c(A) h(X) 
h(X) g(X) 



/r(i;A,ee- 



^ 3 (2A) xH 



/ia(2A) 



/i" ; (l;A,ee- 



(71) 
(72) 

(73) 

(74) 
(75) 

(76) 
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Also one expects that the operators .612(A), .813(A) and .823(A) play the role of creation 
operators over the reference state |0i). Therefore it is natural to seek for other eigenvectors of 
ti(A) as linear combinations of products of these creation fields acting on |0i). This is done by 
exploring the commutation rules between the diagonal Ai(X) and the creation fields which can 
be derived from the boundary Yang-Baxter algebra (JUJ. A careful analysis of these relations 
reveals us that the construction of the eigenvectors can be based on either .812(A) and -813(A) 
or 523(A) and .813(A) pair of fields rather than on arbitrary combination of the three possible 
creation operators. We remark that this redundancy is not particular of this system, but it 
is a general feature of the algebraic Bethe ansatz framework developed in ref . fI5\ l2l>] for a 
large family of integrable vertex models with periodic boundary. This formalism has been 
generalized by Li et al ref. [2? to include vertex models with open boundaries based on ideas 
first envisaged by Fan [21] and later extended for systems solvable by nested Bethe ansatz 
|28j . We also recall that in the context of three state vertex models this approach was recently 
reviewed in ref. [29 . Considering that such algebraic framework has already been described in 
these references, we shall not repeat the details here, and in what follows we will present only 
the main results for the eigenvectors and the eigenvalues. Here we consider that the eigenvectors 
will be constructed in terms of a linear combination of products of the creation fields .812(A) 
and .813(A) acting on the vector |0i). It turns out that the eigenstates r0 n (Ai, . . . , A n )) form a 
multiparticle structure and they can be constructed as 



such that the vector y? n (Ai, . . . , A n ) satisfy a second order recursion relation of the form 



"0n(A 



. . . , A n )) = V ? n(Al, • • • , A n ) |0i) 



(77) 



y?n(A 



. . . , A n ) 



-8i2(Ai)y? n _i(A 2 , • • • , A n ) 



n 



+ -8i3(Ai) V 9 «-2(A2, • • • , Aj_i, Aj + i, . . . , A n ) 
i=2 



x (if (A 1; . . . , \ n )Ai(Xi) + lf(Ai, . . . , K)A 2 {xSj , (78) 



Here we are assuming the identification |?/>o) = |0i). The functions 1^ (Ai, . . . , A n ) and 
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r^(Ai,...,A n ) are given by 



r^(Ai, . . . , a„) 

and 

ri* (Ai, . . . , A n 



d(Xi - A^ _ 



i-l 



6(Ai + a,; 



p(Ai,AOn 



h^{Xj — A,;) 



n 



i-l 

n 



q(A fc - Xj)b(X k + Xj) 
2 a(Aj - Ai)e(Aj - Ai) ^1 l(\ k _ Ai)a(A fc + Ai) ' 

(79) 
^(A fc + A^ 



h^{Xj — Aj) 



n 



/i 4 (Aj — Afc) 



rf(Ai + A^ 

6(Ai + Aj) A = -^ a(Xj - Aj)e(Aj - Aj) ^ 6(Aj - A fc )e(Aj - A fc ) a{X k + Aj)6(A fc + Aj) 



50) 



With p( X ,y) = §^±4 - ^M|^ + 4 and /, 4 ( A ) 



MA) d(A) 
J(A) e(A) 



e(x - y) a(2y) d(x - y) 
The action of the doubled transfer matrix t\(X) on the state |V>(Ai, . . . , A n )) is performed 
relying on similar data for the (n — 1) and (n — 2) particle states and with help of mathematical 
induction. Adapting the discussion of refs. fZT\ l2Hj to our case we can infer that the eigenvalue 
expression is 



MA) 



= 4 +) (X)coi-\x) 
+ J 2 +) (X)J-\X) 
+ 4 + \X)J 3 ~\X) 



q(A) 2 
Ci(A) 

RA) 2 
Ci(A) 

e(A) 2 
Ci(A) 



L n _ 



a( Aj . - A)6(A 3 - + A) 
J- 6(Aj - A)a(Aj + A) 

e(A - X 3 )b{X - Xj)a(X + A i )6(A + A,) 
6(A - A i )/i5(A + Aj) 



n 

L n 

n 



L n 



n 



Le(A-Aj)e(A + Aj)&(A + Aj) 



51) 



where /2.5(A) = 
ansatz equations 



6(A) d(A) 
J(A) 6(A) 



and provided that the rapidities Aj satisfy the following Bethe 



«(Aj) 

5(Ay) 



2 M +) (A>^(Aj) [&(2Aj)] ; 



n 



6(Aj - Xi)h 2 (Xj + Aj 



c4 +) (Aj>W(Aj) /i 2 (2Aj) g(Aj - A,,) [6(Aj + A,)] 



2 • 



52) 



Now we are almost ready to get standard expressions for the eigenvalues and Bethe ansatz 
equations. By introducing the new set of variables Aj = Ai — 77 and performing many simplifi- 
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cations in the functions entering Eqs. f)8H l%2"|) we conclude that the eigenvalues are 



Ai(A) 



(2A + 3r))(\ + eitfL - f)(A + er,C- + jf)(A - i£ + - §)(A - ^ 



(A + 2r]f 
Ci(A) 



L n 



(2A + r?)?? 4 
(Aj-A + 77) (Aj + A - 77) 



n 



+ 



(A + eqt- - 1)(A - eitf_ + ?)(A - rtf+ - g)(A + rtf+ + f] 



-1 L n 



Ci(A) 



n 



(Xj -X + rj)(X + Xj- rj) (A - Xj + 277) (A + Aj + 277) 
L (A ._ A _ 77 )( A + A . +?7 ) (A - Aj) (A + A,) 



+ 



(2A - V )(X - «£. + 1)(A - ei£_ + f)(A + + g)(A + 77^ + 



A(A-r?) 



Ci(A) 



n 



(2A + r/)77 4 

(A - Xj + 2ri) (A + Xj + 2 V ) 
(A - Xj) (X + Xj) 



while the Bethe ansatz equations for the new rapidities Aj becomes 

'A, + v £+ + *\ W (Aj - Aj + 77) (Aj + A, + 77) 



"Aj +77" 


2L 


/ Xj - ei\i- + 


.Aj - 77. 




VAj- +£T7£_ - 



Aj 



n 



2 / if ( A j _ Ai - 77) (Aj- + Xi - 77) ' 



13) 



(84) 



where £± are taken from Eq.(|21|) with S = 1. 

We finally remark that the results of this subsection offers us in principle the basis to solve 
the 0(3) non-linear sigma model with non-diagonal open boundaries. Due to the isomorphism 
0(3) ~ SU (2)2 the elements of the operator (JHHJl can indeed be interpreted as the scattering 
amplitudes of the S-matrix associated to the 0(3) field theory [3*0] . One expects that similar 
relation is also valid for the boundary scattering matrices [HI]. In this case, we need to adapt 
our results to include the solution of the eigenspectrum of an open transfer matrix in the 
presence of inhomogeneities, following for example the lines of ref.|32j. It would be interesting 
to exploit this possibility to determine the effects of the boundaries in the physics of the 0(3) 
model. 
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4.3 The spin-,? solution 



The classical analogue of the solvable spin-S* XXX model is the 2S + 1-state vertex model 

© having the total number of §(2S + 1) [2(2S + l) 2 + 1] non-null Boltzmann weights. The 

■(- 



transformed upper triangular Kg \X) matrix corresponding to the left boundary in tg{X) is 



K ( S -\X) 








\ 







ft\S;X,e^ 



\ 











ftsUS;X,e^) J 
(85) 



where * denotes non- vanishing values that can be directly determined from Eq. (j31|) . 

To implement the quantum inverse scattering framework we will represent the doubled 
monodromy matrix by the following structure 



/ 



rf(x) 



A 1 (X) 
C 21 (A) 



B l2 (X) 
A 2 (X) 



Bips+i) (A) 
(A) 



\ 



B- 



2(25+1) 



B 



2S(2S+1) 



(A) 



16) 



\ C(25+l)l(A) C'(25+1)2(A) ■ • ■ C'(25+l)2s(A) ^2S+l(A) J 

The next step in the algebraic formulation consists in determining the action of the 7"^'(X) 



elements on the reference state \0s) which helps us to distinguish creation and annihilation fields 
as well as to reformulate the eigenvalues problem in terms of appropriate linear combinations 
of diagonal fields. To perform that we need to know certain commutation relations between 



the operators T^\X) and T^^A) . This can be obtained by using Eq.(j2J) with fi = —X 
[21] to get the following general matrix relation 

Tf\-X)Y 1 £ ( £ , (2X)T} S \X)=TI S \X)£[*\2X) [t 2 {s \-X) 

By applying both sides of relation (jHTjl on the pseudovacuum \0g) and by taking into account 
the upper triangular property of both C^(2X) and K s \X) when acting on the state |0s) we 



-i 



87) 
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conclude that all the fields C a p{X) are annihilators, 

C a fs(X) \0 S ) = (88) 

while B a p(X) acts as creation fields upon |0g). 

We also see that Ai(\) |0g) for i = 2, . . . , 2S + 1 turns out to be proportional to many 
distinct bulk terms of the form [tj(A)] 2L since it involves the action of upper elements of the 



operator 
are 



(S), 



-A) 



on |0s). In the specific case of the XXX — S model the functions U(X) 

A + r/k — rjS 



U(X) = (X + 2 V S) ]"[ 



A + rjk + rjS ' 



9) 



k=S-i+2 

As remarked in previous sections this difficulty can be circumvented by writing the fields 
Ai(X) as linear combinations of new operators Aj(A) such that their action on |0g) is propor- 
tional only to [t«(A)] 2L term. The solution of this problem involves a considerable amount of 
algebraic work but the final answer can fortunately be given in terms of the determinants of 
certain j xj matrices that shall be denoted by M^\' (A). Its elements are determined in terms of 



25+1 



the entries of the C^(X) operator. More precisely, by writing C^(X) = -R^(A)e c b (g) e, 
we find that such linear combination is 



ad 



abcd=l 



A(X) = 

3=1 



M«(2A) 



M«(2A) 



MX) 



(90) 



where the j x j matrix M-^(A) is given by 



( Rlj(X) Rlj(X) 



Mff(A) 



^2>(A) Ri]i{X) 



2,2, 



R 



■3-1,3-1 
2,2 



(A) RUX) 



(91) 



By using the relation (j9~Uj) and the action of all Aj(A) on the reference state we find that 



A(X) |0s> = /^^(A) 



r ^ 2 (A) 
Cs(A) 



-i L 



|0s) 



(92) 
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where 



^ H (A) = (e) 



25 



i K-1, fc (2A) 



/<->(£; A, e£_) 



while the entries of a second j x j auxiliary matrix (A) are given by 



/ pi 



m[- } (a) 



fl},i(A) i#?(A) 



l£i M-1 (A) i?iV J+i (A) ^Y' l+1 (A) 



#2,2 (A) #2^) 

V^(A) J#(A) 



-•2,2/ 



1-1,1-1/ 
2,2 



(A) Ri^ j+ \X) R 



i+l,t+l 
2,2 



i— 1 , 



(A) /<;'' ir, (A) a 



i+l,t+l 



(A) 



(A) 



^i(A) 

#2,2 (A) 

^i(A) y 



(93) 



\ 



(94) 



Equipped with Eq.(|9(J|) one can now write the eigenvalue problem (127)) in terms of the new 
diagonal fields A (A), namely 

25+1 

\ , . I A I 

(95) 



J2 ^\X)A k (X) |0) = |v) 



k=i 



where 



(+) 



25+1 



i=k 



Af«(2A) 



MiJ(2A) 



P5 +) 



(96) 



At this stage we would like to emphasize that the role construction presented above is 
applicable to any multistate vertex model whose Boltzmann weights are invariant by one U(l) 
charge conservation symmetry. In order to get manageable expressions for the eigenvalues, 
however, one still needs to carry on cumbersome simplifications on the general formulae given 
in Eqs. (j93l96|) . In the case of the XXX — S model we are able to show that all the contributions 
to J\ (A) miraculously factorized in the following product forms 

w i (±) (A) = r 4 (±) (A)xi ±) (A,e±) 



where 



.(+) 



(A) = n 



k=l 



2A + [25 + 3 - i - k] rj 
2A + [2 - k] 77 ; 



7i<">(A) = (ef s J] 



25 TT 2X + [2 - 2i + k}r] 



k=i 



2A + [1 + k - i] rj ' 



(97) 

(98) 
(99) 
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and 



x! +) (A,e + ) 



25+1- 

n 

3=1 
25 

n 



? n 1 ■ A 
£+ + S+-~J-- 

2 ?7 



i-1 



' ? n 1 ■ A 

e£_ + 5 + - - j + - 
2 ?7 



n 

J j=l L 

25 

n 

J j=25+2-i 



F a 3 ■ A ' 
£+ + S+--J + - 
2 77 



1 



(100) 
(101) 



Before proceeding with further results we stress that the above explicit expressions for 
u\ (A) with arbitrary S are novel in the literature since they were unknown even in the 
case of diagonal boundaries Now we reached a point in which we have gathered the 

basic ingredients to start an algebraic Bethe ansatz analysis of the eigenspectrum of £~s(A). In 
particular the vector |0g) is itself an eigenstate of £~s(A) with the eigenvalue 



25+1 



(102) 



The other eigenvectors of £5 (A) are looked as states created by the action of the fields 
B a /3(X) on the reference state |0g). A single particle excitation is made by lowering the value 
of the azimuthal spin component by an unity on the ferromagnetic pseudovacuum |0g). From 
the point of view of the algebraic Bethe ansatz framework this excitation can be represented 
by Bjj + i(Xi) \0s) for any j = 1, . . . , 2S. As far as commutation relations are concerned we find 
that it is simpler to choose the one-particle state as 



|^i(Ai)> = B 12 (Xx) \0 S ) 



(103) 



The action of the double transfer matrix £5 (A) on this state can be computed with the aid 
of the commutation relations between the fields A* (A) and .812(A) that can be obtained from 
the boundary Yang-Baxter algebra (J41j) . In Appendix B we present details of our analysis of 
the one-particle eigenvalue problem for S = §. This study together with the previous results 



for S = 1 [271 129j and the help of mathematical induction lead us to the following general 
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expression 



25+1 



l&(*o> = E 



i=i 

2S 



Cs(A) 



^ +) (AVr J (A)g i (A,A 1 )|^ 1 (A 1 )> 



,(-) 



+ ^^ +1 (A) ? P(A,A 1 )A 1 (A 1 ) + gf ) (A,A 1 )A 2 (A 1 ) |0 5 > , (104) 



i=l 



where function Qi(A, A-,) is given by 

' R i,l(^j - A)-R?; 2 (Aj + A) 

5 ^ = 



Qi{\ Aj) 



"(A -A,) 



-RiLi(A + Aj 



R. 



^2,i(A + 



2,t 



i^}(A + A,)^(A + A 



it 



^(a - Aj).r5+3:(a - Aj) 

^i,2s(A + Aj) i^ 2S 
^2,25+1 (A + Aj) #2,25+1 (A + Aj 



8,1, 



2,... 



25+1,2 
2,25 



i(A - A, 



25+1,2 (a + A i ; 



, ^i,25+i(A - Aj,) R(%(^ + Aj)-Ri^5 + ^(A + A 



25,1, 



,25+1,1, 



25 + 1 



(105) 

From ()104|) we see that the unwanted terms proportional to Bu + i(\) can be eliminated by 
imposing that the rapidity Ai satisfies the following one-particle Bethe ansatz equation, 

^(AOfM^Ai) ^(A.Ai) 



* 2 (Ai 



1 2S. 



(106) 



We note that though the expressions for g} 1J (A, Ai) and qf\\ Ai) have in general a very 
involved dependence on the i-th index, see for instance Appendix B, we have found out that the 



ratio 



gj 2) (A,Ai) 



does not depend of such index. Its expression for arbitrary S, coming directly 



from the commutation rules, involves many complicated terms and it has been collected in 
Appendix D. It turns out, however, that it is possible to carry out further simplifications in 
equations ()D.l|D.2j) thanks to several identities between the Boltzmann weights R^ b (X). This 



J 2 ) 



(A,Ai) 



also leads us to conclude that the ratio %n 



does not depend on the spectral parameter A. 



This is consistent to what one would expect from a standard Bethe ansatz analysis and the 
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simplified expression for such ratio reads 



^ +) (Ai) 



e(Ai 



(107) 



g^CA.AO cv[ +) (X 1 ) 
where for later convenience we define function 0(A) separately, namely 

Ki(A)] 2 

Putting all these results together we find that |i/>i(Ai)) is an eigenvector of ts(X) with 
eigenvalue Ax(A, Aj.) given by 



A«(A,A X ) 



-i L 



t=l 



.Cs(A) 



^ (+) (A)^- ) (A)Q J (A,A 1 ) 



(109) 



provided that the variable Ai satisfies the restriction 

.t 2 (Ai 



e(Ax). 



(110) 



wJ +) (Ai)^ _) (Ai) 

Here we remark that the equation (jllOj) is equivalent to the condition of analyticity of 
Ai(A, Ai) as a function of the rapidity Ai. This fact is indeed an extra verification of the 
validity of our Bethe ansatz analysis. 

We now turn to the analysis of the two-particle state. In this case one expects that this 
state should be given in terms of two linearly independent vectors -Bi2(Ai)-E>i2(A2) \0s) and 
-Bi3(Ai) |0s). Previous experience in determining two-particle states [2U 123 I2H] suggests us 
to look first for the commutation rule between the fields -Bi 2 (Ai) and -B 12 (A 2 ). To avoid 
overcrowding this section with more heavier formulae we have exhibited this relation for S > 1 
in Appendix D. From equations (j9*Uj) and the observations made in Appendix D we clearly see 
that the state 

-^i! 2 (Ai + A 2 ) 



|$ 2 (Ai,A 2 )) = J B 12 (A 1 ) J B 12 (A 2 ) + S 13 (A! 



+ 



^i' 2 (Ai + A 2 ) 


M$(2\ 2 ) 


rI^(Xi 


^i' 2 (Ai + A 2 ) 


MW(2A 2 ) 


#i>(Ai 



^i,' 2 (Ai + A 2 ) 

\ \ n3,l/i 



>2,1, 




|0s>, (HI) 
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is symmetric under the exchange of the rapidities Ai and A2. In other words we have that 

|^ 2 (A 1; A 2 )) = Z S (X U A 2 ) |Vi 2 (A 2 , AO) (112) 
where Zs(\i, A 2 ) is the following function, 



Zs(Xi, A 2 ) 





i?l' 2 (A - 


AO 


i^J(A-A0 




-R 22 (A - 


AO 


R l',3 (A - AO 



(113) 



^(A + AO ^(A-AO^kA-AO 
This state is therefore an educated ansatz for the two-particle vector for general S > 1. 
Note that it reproduces the previous state for S = 1 [2H 123 I2H] and in Appendix B we 
have presented all the needed evidences that it is indeed a suitable eigenvector for S = §. 
The corresponding eigenvalue can be calculated by keeping only the terms proportional to the 
vector i?i 2 (Ai)i?i 2 (A 2 ) coming from the first part of the commutation relations between the 
fields Ai(X) and _Bi 2 (Aj)- Taking into account our previous experience with the one-particle 
state and the structure of the commutation rules discussed in Appendices B and D we find 
that 



A?'(A,{Ai,A 2 }) 



25+1 

£ 

i=l 



Cs(A) 



n=2 



^ (+) (AK ( - ) (A)n^(A,A J 



;ii4) 



The associated Bethe ansatz equations are expected to be the condition on the rapidities 
such that the residues at the simple poles A = Ai,A 2 present in functions Qi(X, Xj) vanish. 
This condition is equivalent to the following system of equations 



j / 



2L , ,(+) 



c^(A>j _) (A0 



71=2 



j 



n 



(115) 



By the some token, one expects that general multiparticle states can in principle be 
constructed in terms of a recurrence relation of order 2S that involves the creation fields 
-^ij(A) j — 2, . . . , 2S + 1. The precise structure of such relation for arbitrary 5* has however 
eluded us so far. This by no means prevents us to propose general expressions for the corre- 
sponding eigenvalues and Bethe ansatz equations. In any factorizable theory, it is believed that 
the two-particle results already contain the main flavour about the content of the spectrum. 
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This means that the expressions (| 114)1 and pi 5)1 are expected to be valid for general values 
of n < 2LS. Considering these observations and after working out the explicit expressions for 
functions Qi(X, Xj) we find that the n-particle eigenvalue A§ (A, {A;}) is given by 

A?(A,{A,}) 



n +] l r ) " ^ 

Ps Ps i=l 



LCs(A)J 

[A - Xj + rj(S + 1)][A - A, - rjS] [A + A, + r)(S + 1)][A + A, - rjS) 

1 [A -A 



n 



X " J - A, + v(S + 2 - i)][A - Xj + V (S+1- i)} [A + A, + ij(S + 2 - i)][X + Xj + V {S + 1- i)} 



(116) 



while the Bethe ansatz equations are given by 



Xj + 7]S 


2L 


(Xj 


- + 


Xj — T]S 




\Xj 





I \, + Vt+ + !\, A (Aj - Ai + 77) (A, + A, + 77) 



>/ , (A, A, '/) (A/ • A, ,/)• 1 '' 



where we have performed the displacement Aj —>■ Aj — i]S in order to bring the Bethe ansatz 
equations in a more symmetrical form. 

At this point it should be emphasized that the right hand side of the Bethe ansatz equations 
pi7)) depend on both the spin S and the off-diagonal elements c±, d± through the renormalized 
variable £± defined in Eq. (|21j) . We also mention that we have verified numerically for several 
values of L and S that the equations pi6H17)) indeed reproduces the ground state and few 
low-lying excitations of the double- row transfer matrix t$(X). In particular, we have been able 
to check the completeness of the Bethe ansatz solution for L — 2 up to S — |. Finally, we 
remark that the final results for the eigenvalues (J114)) and Bethe ansatz equations pi 5)) are 
expected to be valid for any integrable vertex model whose underlying .R-matrix possesses an 
unique U(l) charge symmetry and the invariance (J5I6)) . 



5 Conclusions 

The purpose of this paper was to solve the integrable XXX-S 1 Heisenberg model with open 
boundary conditions by means of the quantum inverse scattering approach. We first argued 
that the corresponding fT-matrices are diagonalizable by special similarity transformations 
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without a dependence on the spectral parameter. This fact together with the property of 
reversing gauge transformed Boltzmann weights leads us to an eigenvalue problem with only 
one non-diagonal effective i^-matrix. In the cases when such i^-matrix are either upper or 
lower triangular we managed to present explicit expressions for the eigenvalues of the doubled 
transfer matrix operator £s(A) as well as the associated Bethe ansatz equations for arbitrary 
values of the spin-S'. This condition was shown to be equivalent to two possible constraints 
between the four off-diagonal boundary parameters, leading us with five free parameters out 
of six possible ones. 

We hope that the ideas developed in this paper will be also suitable to solve a broad class of 
isotropic integrable systems with non-diagonal open boundaries. In fact, the method devised 
here has been first applied to the fundamental SU(N) isotropic vertex model under more 
restrictive open boundary conditions [33]. We expect that the nested Bethe ansatz approach 
could be further generalized to tackle effective triangular .fT-matrices which will provide us the 
solution of the associated doubled transfer matrix operator with fewer constrained boundary 
parameters as compared to that presented in ref.[33]. We also hope that other vertex models 
based on higher rank symmetries such as O(N) and sp(2N) Lie algebras could be dealt by 
the framework discussed in this work. This assumes that certain classes of non-diagonal K- 
matrices of these vertex models can be classified in terms of similarity transformations that 
are itself symmetries of the corresponding ^-operator, acting on spectral dependent diagonal 
solutions for the reflection equation. This would means that our observation of section (J2J) for 
SU(2) could be generalized to other Lie algebras as well. We plan to investigate such rather 
interesting possibility in a future work. 
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Appendix A: The i^-matrix properties 



In this Appendix we briefly summarize the explicit expressions of the if-matrix elements 
satisfying the reflection equation (JHJ) for S = 1 and |. For S = 1 [TT| the corresponding matrix 
is given by 



( MA) MA) MA) ^ 



^i(A) 



fe 21 (A) k 22 {\) k 23 (X) 

y fe 3 i(A) /c 32 (A) a; 33 (a) y 

where the elements fcy(A) are given by 



feu(A) 
MA) 
MA) 
fei(A) 
MA) 
MA) 
MA) 
MA) 
MA) 



2 77 



1 

4 

c 



„ 1 A\ cd 
2£+- + - +- 



»7 



2^ 

~2 



1 A\ A 

2 T] J 7] 

X\ X 



d 



2^2 
1 

~4 

c 

2v/2 

£ 

4 

d 

2V2 
1 

~4 



2£ 
2^ 

1 

2 _ 



77/77 
1 



A\ A 
2 77/77 
_ A 
77 



2£ 
A\ A 



+ 



cci 



2£4 
2^ 



1 

2 77/77 

77/ 77 

1 A\ A 

2 77/77 

1 _ A 

2 ~ 77 



2 77 



+ 



erf 



On the other hand for S = | we have 



i^3(A) 



A;n(A) fc 12 (A) MA) fcu(A) 



MA) MA) MA) MA) 

fe 31 (A) MA) MA) MA) 

^ MA) MA) MA) MA) J 
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(A.l) 



(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 
(A.9) 
(A.10) 



(A.ll) 



where the elements fcy(A) are given by 



MA) = £ 



V 



77 



+ ^(3^+-) (3^-1 + 



fci 2 (A) = 
MA) = 
MA) = 
MA) = 

Ma) = 
+ 

MA) = 
hi (A) = 

Ma) = 

MA) = 
MA) = 



18^ 



erf 1 - 



(\ 



V!U-- ),3e - 1 + - ) -' 



2 ( 3£ - 1 + 
A A A 



A 



3{ + 



A N 



A 



c3/l_AA / AAA 
27 V 2 v) V 



d 



18^ 
2cd 



27 



cd f 1 - - 

3 



+ 3£- 



2(3^ 

A 
77 



77 



77 



A 
77 



+ 



1 

4 + 4 



3f \. A 

77 



2 



77 



3e 



77 



c 

'54 



cd 1- 



1 + 



2A 
77 



c 2 (\ 



A 



9^ V 2 ?7 
d 2 ( 1 A 

9^ V 2 77 



+ 4 



A A A 



1 - 



- (30 2 



3f + 1 - 

77/ 77 

A A A 

3£ - 1 + - - 
77/ 77 



"54 

2cd 
~27 



cd 1 - 



A 



- 3£ + 



+ 



27 V 17 



2A 
77 

4) G 
A 



4 



1 - 



" (30 2 



+ 



3£ 



3£ + 1 - - 
77 

k u (X) = \cd ( 1 - -) - 2 f 3£ - - 

18^3 [ I f|j V V 



i + x - 

77 



4 j 77 2 



A 

■ 

77 



3£ + l- 



77 



77/77 
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MA) 



18-s/3 



A 
r/' 



(A.28) 



18 



77 



1-* 

77 



4M 

27 V V 



77 



3^ + 1 



Next we list the dependence of the off-diagonal coefficients of the transformed 
Kg~\\) on the parameters c± and d±. For S — \ we find 

-e+ [c+cL + c_d + — 2c + d + ] + (c_d + — gLc + ) ^/l + c+d + 



C12 : 

while for S — 1 we have 

(2 + C-d + + c + gL 



«12 



«13 



«23 = 



32v^(l + c+rf+) 
-e + 0° + 24 (d+ 



(-e+2(c+ - c_)v/l + c + rf + - A) , 



d. 



8\/2^1 + c + d 
where A and 0* are given by 



32(1 + c+d+){2 + c+d+ + 2e +v /l + c + d + ) 2 
(2(c + - c-)y/\ + c + d + + c+ a) , 



(A.29) 
i^-matrix 

(A.30) 



(A.31) 

(A.32) 
(A.33) 



A 


= 2c_ — c + (2 — c_d + + c + d), 


(A.34) 


6° 


= c+(d + - d-)y/l + c + d + [cL(2 + c + d+) - d+(2 + c_d + )] , 


(A.35) 


e 1 


= 4c + d + + c_rf + (4 + c_d + ), 


(A.36) 


e 2 


= 2c_d + + c + d + (6 + c_d+), 


(A.37) 


e 3 


= c+d + (8 + c + d + ). 


(A.38) 



Appendix B: One and Two particle states for S = | 

The purpose of this Appendix is to present some of the technical details entering the analysis 
of the one and two particle states for S — §. In order to do that it is convenient to work with 
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a new matrix -R„f (A) = P a bCf b (X) where P ab is the permutator. This matrix plays a direct 

role in the quantum inverse scattering method and Eq.(|41j) is rewritten in terms of R^ (A) as 

1 11 1 

R[ S 2 \u - v) Tf (u)R[ S 2 \u + v) Tf (v) =Tf (v)R[ S 2 \u + v) (u)R[ S 2 \u - v), (B.l) 

- (-) 

In the specific case of a 44 vertex model, the R a l (A) operator can be expressed in terms 
of the following matrix: 



a(A) 


















































6(A) 








e(A) 









































c(A) 








/(A) 








MA) 
































d(X) 








9(A) 








i(A) 








i(A) 














e(A) 








6(A) 









































r(A) 








1(A) 








r(A) 
































91 (A) 








m(A) 








9(A) 








n(A) 
































n(A) 








n(A) 








MA) 














h(X) 








/(A) 








c(A) 
































n(A) 








9(A) 








m(A) 








9i(A) 
































/i(A) 








MA) 








/i(A) 









































MA) 








ei(A) 














iW 








i(A) 








9(A) 








d(X) 
































hi (A) 








nW 








n(A) 









































ei (A) 








MA) 


















































ai(A) 



(B.2) 



In order to solve the one-particle problem one first needs to obtain the appropriate com- 
mutation rules between the fields Ai(u) and E>i 2 (v) coming from the boundary Yang-Baxter 
equation (jB.lj) . Using the symbol [i,j] to represent the i-th row and the j-th column of Eq.(jB~T| 
we conclude that such suitable commutation rules are derivate from the entries [1,2], [2,3], [3,4], 
[2,6], [3,7] and [4,8]. Further progress are made replacing the fields A^u) by A^u) in these 
equations with the help of the relations (J90|) . After several algebraic manipulations we obtain 
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the following structure 

Ai(u)B 12 (v) = a\(u,v)B 12 (v)A 1 (u) + al(u,v)B 12 (u)A 1 (v) + al(u,v)B 12 (u)A 2 (v) 

+ a\(u,v)B 13 (v)C 21 (u) + al(u,v)B 13 (u)C 21 (v) + a\(u,v)B 13 (u)C 32 (v) 
+ a\(u,v)B 14 (v)C 31 (u) + al(u,v)B u (u)C 31 (v) + al(u,v)B 14 (u)C 42 (v) 

(B.3) 

A 2 (u)B 12 (v) = al(u, v)B 12 (v)A 2 (u) + a\(u, v)B 12 (u)A 1 (v) + a 2 3 (u, v)B 12 (u)A 2 (v) 
+ al(u,v)B 23 (u)A 1 (v) + al(u,v)B 23 (u)A 2 (v) + al(u,v)B 13 (v)C 21 (u) 
+ a 2 (u, v)B 13 (v)C 32 (u) + a 2 8 (u, v)B 13 (u)C 21 (v) + a%(u, v)B 13 (u)C 32 (v) 
+ a\ Q (u,v)B 24 (u)C 21 (v) + a 2 11 (u,v)B 24 (u)C 32 (v) + a 2 12 (u, v)B u (v)C 31 (u) 
+ a\ 3 (u,v)B 14 (v)C 42 (u) + al 4 (u,v)B 14 (u)C 31 (v) + a 2 15 (u, v)B 14 (u)C 42 (v ) 

(B.4) 

A 3 (u)B 12 (v) = a\{u, v)B 12 (v)A 3 (u) + a 3 (u, v)B 12 (u)A 1 (v) + a\(u, v)B 12 (u)A 2 (v) 
+ a 3 4 (u, v)B 23 (u)A 1 (v) + a 3 (u, v)B 23 (u)A 2 (v) + a 3 (u, v)B 34 (u)A 1 (v) 
+ a 3 7 (u,v)B u (u)A 2 (v) + a 3 8 (u,v)B 13 (v)C 21 (u) + a 3 9 (u,v)B 13 (v)C 32 (u) 
+ a 3 w (u,v)B 13 (v)C 43 (u) + a 3 u (u,v)B 13 (u)C 21 (v) + a\ 2 {u,v)B 13 {u)C 32 {v) 
+ al 3 (u,v)B 24 (u)C 21 (v) + a 3 u (u,v)B 24 (u)C 32 (v) + a? 5 («, v)B u (v)C 31 (u) 
+ al 6 (u,v)B u (v)C 42 (u) + al 7 (u,v)B 14 (u)C 31 (v) + a? 8 (u, v)B 14 (u)C 42 (v) 

(B.5) 
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A 4 (u)Bi 2 {v) = a\(u, v)Bi 2 {v)A A (u) + a 2 (u, v)B 12 (u)A 1 (v) + a^(u, v)B 12 (u)A 2 (v) 
+ a\(u, v)B 23 (u)A 1 (v) + 0,5(11, v)B 23 (u)A 2 (v) + a^u, v)B u (u)A 1 (v) 
+ aj(u,v)B u (u)A 2 (v) + al(u,v)B 13 (v)C 21 (u) + ag(u, v)B 13 {v)C 32 (u) 
+ a^ (u,v)B 13 (v)C i3 (u) + a\ x (u,v)B 13 (u)C 21 (v) + a\ 2 (u, v)B 13 (u)C 32 (v) 
+ af 3 (u, v)B 2i (u)C 2l {v) + af 4 (u,v)B M (u)C 32 (v) + af 5 (u,v)B u (v)C 31 (u) 
+ a* 6 (u, v)B u (v)C A2 (u) + a\ 7 (u, v)B 1A (u)C 3 x(v) + af 8 {u, v)B 14 (u)C 42 (v) 

(B.6) 

Before proceeding we would like to remark that several identities between the Boltzmann 
weights have been used in order to obtain relations f)B.3IIB.6|) . We also note that many of the 
coefficients a\[u, v) are proportional to annihilation operators and not all of them are relevant 
in the calculations. In Appendix C we have listed only those that indeed play an important role 
in our analysis since in general they are sufficiently cumbersome. By applying Eqs. (jB.3HB.6j) on 
the pseudovacuum O3 ) we see that one can rearrange the action of the double transfer matrix 



£3 (A) on the one-particle state Bi 2 (\i) 



O3 ) as in Eq. ()104|) . Furthermore, it turns out that the 



functions Qi(X, Ai), ^(A, Ai) and qf\\ \\) can therefore be explicitly read off, namely 

Qi{\ Ai) = oi(A,A0 (B.7) 
4 

qP(\ Ai) = ^^(A^AA) (B.8) 



3=% 



q?\\\i) = J2^ + \X)4 i+1 (X,X 1 ) (B.9) 



3=1 



where % = 1, . . . , 4 and function t<4 (A) has been defined in Eq. (}96|) . 

(2),, , \ 

As mentioned in the main text the ratio g ' m ' is independent of the i-th index and of 
the spectral parameter A. In our case this ratio is given by 

of\\M 4 +) (AQ / a(Ai)/(Ai) - b(Ai)b(Ai) \ 
g «(A,A0 w «(A0V [e(Xi)} 2 J ( ' Uj 

Next we turn to the two-particle state. In order to obtain an ansatz to this vector we have 
considered the commutation rules [1,3] and [1,6] coming from Eq. (jB.l|) . Acting these relations 
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on 



O3 ) leads us to the following expression 



B l2 (u)B l2 (v) + B 13 (u) (a 2 (u,v)A 2 (v) + a i (u,v)A 1 (v] 



■1 



Zs (u, v) 

2 v ' 



B 12 (v)B 12 (u) + B 13 (v) (a 3 (u, v)A 2 (u) + 014(11, v)A 1 (u) 



where functions ai(u, v) for i — 1, ... ,4 and Zz (u, v) are given by 

f(u + v)b(2v) f(u — v)h(u + v) 



a>i{u } v) 



e{u + v) a(2v) h{u — v)e(u + v ) 



(B.11) 



(B.12) 



f(u + v) f(u + v) 

a 2 (u,v) = — — ■ — -, a 3 {u,v) 



e(u + v ) 



e(u + v) 



e(u + v) a(2u) \ h(u — v)l(u — v) — f(u — v)r(u — v ) J e(u + v) 



h{u — v)l(u — v ) — f(u — v)r(u — v) \ e(u + v) 



Zziu, v ) — , 

2 ' \ a{u — v)h{u — v) J e(u + v) 

From Eq. ()B.ll|) it follows that an appropriate two-particle state should be 

<MAi, A 2 )> = fe 12 (A 1 )S 12 (A 2 ) + S 13 (Ai) (a 2 (Ai, A 2 )1 2 (A 2 ) + ^(Ax, A 2 )2 1 (A 2 ) 



(B.13) 



(B.15) 



(B.16) 

since it is symmetric |^ 2 (Ai, A 2 )) = Zz(\i, A 2 ) |^i(A 2 , Ai)) under the exchange of rapidities. 

The next step is to solve the eigenvalue problem for the two-particle state (jB.16|) . In order 
to do that we need extra commutations rules between the fields Ai(u) and Bi 3 (v), Bi 2 {u) 
and Bjj + i(v), Cj + ij{u) and B 12 (v). In the case of the fields Cj + ij{u) and Bi 2 (v) the rules 
comes from the entries [2,5], [3,6] and [4,7] of Eq. (jB.l|) and the ones for the other operators 
are obtained from [1,3],[2,4];[1,6],[2,7],[3,8];[2,10],[3,11] and [4,12] entries. After long algebraic 
manipulations we are able to obtain the following expressions 



Ai(u)B l3 (v) 



Oa 



b\(u,v)B 13 (v)A i (u) 



03) + unwanted terms i — 1, 



(B.17) 
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C 21 {u)B 12 {v) 



c l 22 (u, v)A 2 (v)A 2 (u) + c l 21 (u, v)A 2 (v)A 1 {u) + c\ 2 (u, v)A 1 (v)A 2 (u) 



+ c\ 1 (u,v)A l (y)A l (u) 



Os 

2 



(B.18) 



C 32 {u)B l2 {v) 



2 



c 2 23 (u,v)A 2 (v)A 3 (u) + c 2 22 (u,v)A 2 (v)A 2 (u) + c 2 21 (u,v)A 2 (v)A 1 (u) 



+ c 2 13 (u, v)A 1 (v)A 3 (u) + c 2 12 (u, v)A 1 (v)A 2 (u) + c 2 n (u, v)A 1 (v)A 1 (u) 



(B.19) 



C A3 {u)B 12 {v) 



= c 3 24 (u, v)A 2 (v)A A (u) + c 3 23 (u, v)A 2 (v)A 3 (u) + c 3 22 (u, v)A 2 (v)A 2 (u) 
+ cl x {u,v)A 2 {v)A x {u) +4±{u,v)A 1 {v)A 4 {u) + c? 3 (m, w)l 1 (t»)2 3 (M) 



+ «)v4i(v)A 2 (ii) +c? 1 («,v)Ai(v)Ai(«) 



(B.20) 



B 12 (v)B 12 (u) 



B l2 {v)B 23 {u) 



2 



ci^w, i')-^i3( t ') y 42(w) + v)B V3 {y)Ai{u) + unwanted terms 



0| 



(B.21) 



d 2 3 {u, v)B 13 (v)A 3 (u) + d 2 (u, v)B l3 {v)A 2 {u) + u)Si 3 (t;)Ai(u) 



+ unwanted terms] 



0: 



(B.22) 



B 12 (v)B u (u) 



(dl(u, v)B 13 {v)A 4 {u) + d 3 3 {u, v)B 13 {v)A 3 {u) + d 3 2 {u, v)B 13 (v)A 2 (u) 



+ c/^(m, v)B V3 {y)Ai{u) + unwanted terms 



a 

2 



(B.23) 



where by "unwanted terms" we mean those that do not give contributions proportional to 
|^ 2 (Ai,A 2 )). The functions b\(u,v), c^(u,v) and d\{u,v) are once again very involved and 
have been presented in Appendix C. 

We have now the main ingredients to study the action of the operators Ai(X) on the two- 
particle state W 2 (\i, \ 2 )\ Taking into account the commutation rules Eqs. (jB.3|) - (|B.6|) and 
(jB.17j) - ()B.23j) and after some algebra we conclude that the two-particle wanted terms have the 
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following structure 
ts(A) 



Pi 



|&(Ai,A 2 )) = J B 12 (A 1 ) J B 12 (A 2 )^^ (+) (A)(n«UA,A,))A(A) 



Ai=2 



i=l 



U'=l 



+ £x 3 (AiM 4 (A) (Af (A, {A,})1 2 (A 2 ) + A?(\, {A,})ix(A 2 ) 

+ B 13 (X 1 )A 3 (X) (Af (A, {A,})2 2 (A 2 ) + A^(A, {A,})2 1 (A 2 ) 

+ Bx 3 (AiM 2 (A) (A? (A, {A,})2 2 (A 2 ) + A 21 (A, {A i })2 1 (A 2 ) 

+ B 13 (\ 1 )A 1 (X) (A 2 2 (A, {A,})2 2 (A 2 ) + A^(A, {A i })I 1 (A 2 ) 



-unwanted terms 



(B.24) 



where functions A 2 fc (A, {A«}) are given by 



A 2 2 (A, {A,}) = (&?(A, Ai)a 2 (Ai, A 2 ) + <4(A, Ai)4(A, A 2 ) + a*(A, Ax)a 4 7 (A, A 2 )d|(A, Ai) 

+ 4 +) (A) f a? (A, A 1 )4(A, A 2 ) + a?(A, Ai)a*(A, A 2 )^(A, Ax)") (B.25) 



A* 1 (A, {A,}) = 4 +) (A) (^(A, AOa^Ax, A 2 ) + af (A, A^c^A, A 2 ) + a?(A, Ai)ag(A, A 2 )^(A, Ax) 
+ 4 +) (A) f a? (A, A 1 )c? 4 (A, A 2 ) + a?(A, AxK(A, A 2 )t^(A, Ax)") (B.26) 



A 2 2 (A, {A,}) = 4 +) (A) ^ (A, Ax)4(A, A 2 ) + a|(A, A 1 )c 2 3 (A, A 2 ) + a*(A, Ai)a*(A, A 2 )d*(A, Ax) 
+ a*(A, Ai)af (A, A 2 )d 2 (A, Ai)^ + 4 +) (A) ffef (A, A 1 )a 2 (A 1 , A 2 ) + a? (A, Ai)4(A, A 2 ) 



+ ajj(A, Ax)c 2 3(A, A 2 ) + a?(A, Ai)a?(A, A 2 )djj(A, Ax) + a?(A, Ax)a|(A, A 2 )d 2 (A, Ax) 
+ 4 +) (A) (a 2 7 (A, Ax)c 2 3 (A, A 2 ) + a 2 (A, Ax)a 2 (A, A 2 )rf 2 (A, Ax)) (B.27) 
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A 31 (X, {Xi}) = J+\x) \aU\ Ai)4(A, A 2 ) + a$(\, X l )c 2 l3 (X, A 2 ) + a*(A, Xja^X, X 2 )d 3 (X, X l ) 
+ ai(X, X 1 )4(X, X 2 )d 2 3 (X, X,)) + 4 +) (A) (b 3 (X, X 1 )a 1 (X 1 , A 2 ) + a 3 w (X, X 1 )c 3 3 (X, A 2 ) 



+ ajj(A, Ax)4(A, A 2 ) + af (A, Ai)ajj(A, A 2 )^(A, Ax) + a 3 (X, A^a^A, A 2 )d*(A, Ai) 
+ 4 +) (A) f a 2 7 (A, A04(A, A 2 ) + a 2 (A, Ax)a 2 (A, A 2 )d 2 (A, Ai)") (B.28) 



A2 2 (A, {A,}) = 4 +) (A) ( a 4 10 (X, X 1 )c 3 22 (X, A 2 ) + a|(A, A 1 )c| 2 (A, A 2 ) + a|(A, X 1 )c 1 22 (X, A 2 ) 



+ a?(A, AOa^A, A 2 )4(A, Ai) + a*(A, Ai)a£(A, A 2 )d 2 (A, Ai) + a*(A, Ai)a£(A, A 2 )4(A, Ai) 
+ 4 +) (A) ( o? (A, Ai)4j(A, A 2 ) + a 3 9 (X, X 1 )c 2 22 {X, A 2 ) + a|(A, A^c^A, A 2 ) + a?(A, A x ) 



x o 3 7 (A, A 2 )4(A, A x ) + o?(A, Ax)ai(A, A 2 )4(A, A x ) + a?(A, A^a^A, A 2 )4(A, A x ) 
+ 4 +) (A) (fc 2 (A, A 1 )a 2 (A 1 , A 2 ) + a 2 7 (A, A^c^A, A 2 ) + a 2 (A, X 1 )c l 22 (X, A 2 ) 



+ a 2 (A, Ai)a 2 (A, A 2 )4(A, Ai) + a 2 (A, Ai)a 2 (A, A 2 )4(A, Ai) 
+ ^ +) (A) (a\(X, X 1 )c 1 22 (X, A 2 ) + a}(A, Ai)aJ(A, A 2 )4(A, Ai) 



(B.29) 



A 21 (A, {A,}) = u{ +) (X) ( ai (X, X 1 )c 3 2 (X, A 2 ) + a 4 9 (X, X 1 )c 2 12 (X, A 2 ) + a|(A, A 1 )4(A, A 2 ) 



+ a*(A, Ai)ag(A, X 2 )d 3 (X, Ai) + a?(A, A^a^A, A 2 )4(A, Ai) + a?(A, A^a^A, A 2 )4(A, Ai) 
+ 4 +) (A) fa? (A, A!)4(A, A 2 ) + ajj(A, A^c^A, A 2 ) + a|(A, A^c^A, A 2 ) + a?(A, A x ) 



x ajj(A, A 2 )4(A, AO + a?(A, Ai)a|(A, A 2 )4(A, A x ) + a?(A, Ax)4(A, A 2 )4(A, A x ) 
+ 4 +) (A) f& 2 (A, X 1 )a 1 (X 1 , A 2 ) + a 2 (A, A^c^A, A 2 ) + a 2 (A, A 1 )cJ 2 (A, A 2 ) 



+ a 2 (A, A!)a 2 (A, A 2 )4(A, A x ) + a 2 (A, Ai)a 2 (A, A 2 )4(A, A^ 
+ c^ +) (A) faJ(A, AOcl^A, A 2 ) + a\(X, X 1 )a 1 2 (X, A 2 )4(A, A^ 



(B.30) 
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A* 2 (A, = (ai (X, Ai)t| 1 (A, A 2 ) + 4{X, X 1 )c 2 21 (X, A 2 ) + a*(A, A 1 )4 1 (A, A 2 ) 

+ 4{X, Ai)4(A, A 2 )d?(A, Ai) + a}(A, Ai)a|(A, A 2 )d?(A, Ax) + af(A, Ai)ag(A, A 2 )^(A, Ax) 

+ 4 +) (A) f<4(A> Ai)4(A, A 2 ) + a 3 9 (X, Ai)4(A, A 2 ) + ajj(A, Ai)4(A, A 2 ) + a? (A, Ai) 

x 4(A, A 2 )df (A, AO + a?(A, Ai)ajj(A, A 2 )d?(A, A x ) + a?(A, A 1 )a^(A, A 2 )^(A, Ai) 

+ 4 +) (A) fa?(A, Ai)4(A, A 2 ) + a*(A, A^c^A, A 2 ) 

+ a 2 (A, Ai)a§(A, A 2 )d?(A, Ai) + a? (A, Ai)og(A, A 2 )rf}(A, Ai) 

+ wJ +) (A) ffci(A, Ai)a 2 (Ai, A 2 ) + a\(X, Ai)4(A, A 2 ) + a* (A, Ax)4(A, A 2 )^(A, Ai) ) (B.31) 



A" (A, {A,}) = 4 +) (A) U (A, A!)4(A, A 2 ) + a 4 g (X, X 1 )c l n (X, A 2 ) + a*(A, A 1 )cl 1 (A, A 2 

+ 4{X, Xt)4{X, A 2 )d?(A, Ai) + a}(A, Ai)a|(A, A 2 )d?(A, Ax) + a?(A, Ai)4(A, A 2 )d}(A, Ai) 

+ 4 +) (A) K (A, Ai)ci 1 (A, A 2 ) + al(\, A X )4(A, A 2 ) + ajj(A, Ai)4(A, A 2 ) + a?(A, Ax) 

x ajj(A, A 2 )df (A, Ai) + a? (A, Ai)af (A, A 2 )d?(A, A x ) + a? (A, Ai)a^(A, A 2 )rf}(A, Ai) 

+ 4 +) (A) (a 2 (A, Ai)4(A, A 2 ) + a 2 (A, X 1 )c\ 1 (X, A 2 ) 

+ a?(A, Ai)a|(A, A 2 )d?(A, A x ) + a? (A, A 1 )a^(A, A 2 )d}(A, Ai) 

+ ^S +) (A) fftKA, Ai) ai (Ai, A 2 ) + a\(X, Ai)^(A, A 2 ) + aJ(A, A^a^A, A 2 )rf}(A, A x ) ) (B.32) 



It turns out that many identities between the Boltzmann weights can be used in order to 
show the following remarkable property 

n=2 

A l 2 k (X,{X t })=a k (X l ,X 2 )ujl + \X)l[a[(X,X i ) (B.33) 

i=i 

Considering Eq. (jB.24|) . Eq. ()B.33|) and Eq. (j9^j) together it is not difficult to derive the 
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expression 

L /n=2 

^\\)^\\)[\{^{\\ j )\\U^M) 



u{\) _ * 

to W^)) = E 

Pi Pl 1=1 

2 2 



Ca(A) 



i=i 

+ unwanted terms (B.34) 

As a final comment we would like to stress that we have also performed extensive checks 
verifying that in fact the unwanted terms are canceled out provided the rapidities Aj satisfy 
the restriction (jll5j) . 

Appendix C: Auxiliary functions for S = | 

The purpose of this Appendix is to list the expressions of the functions al(u,v), b[(u,v), 
Ci (u, v) and d\{u, v) used in the previous Appendix. To sort notation we shall used the symbol 
u± = u±v and we emphasize that the most complicated functions a^u, v) and af(it, v) have 
been collected at the end of this Appendix: 

a ^ v) = a{u + )e{-u_) (C1) 

^ V) = ~ { a(u + )e(-u.) J (C - 2) 

al(u,v) = -( b -p±) (C.3) 
\a(u + )J 

i / \ a(— U-) f(u+) ,^ 
a\(u,v)= ) Y V + C.4 
a{u + ) e{—u_) 

2 _ (- (6(tx+) 6(tf+)) + a(M+) /(>+)) (%_) l(u-) - f(u-) r(n_)) 

a[u + ) e{U-) e[u + ) a(U-) 

f , U x (b(2u) . ft(u_)6(u + )\ 6(2^ f fe ( 2 ") fe (»+) + fe («-)(^ + ) 2 -^+)^(" + )) \ 

UfeT + e u_ e Lt a (2«)o(«+) + «(«+) e(« + ) J 



.Kfu) "•" e( U _)e( U+ ) i V ' la(2«)o(«+) o(tt+) e(«_) e(«+) J 

w^) = rh — s + 7^~\ 

a(w+J e(— w_J a{2v) 
a(u-)b(u+) (h(u-)l(u-)- f(u-)r(u-)) b(u_) f(u-)r{u + ) 
a(u+) e(w_) 2 h(u-) e(uJ) e(u+) h(uJ) 
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2 b(2u)b(u + ) 6(u_) (6(w+) 2 - a(u+) Z(u+)) 

a 3 m, f) = — , — 7 — H 7 — -— — -— — : (U.7) 

a{2 u) a[u + ) a(u + ) e[U-) e[u + ) 



^«)=-(S) (a9) 



a fi u,iM 



6(«+)6(«_) Z(«_) -/(«_) r(«_)) /(«+) 



a{u + ) e(u_) e(u + ) h(u_) 

:-«-)/(«+) (^ + ^7^)\ | g ( M+ ) (-(/( M _) C ( M _)) + %-)/(«-)) (C10) 



a(u+) e(— w_) / e(w_) e(u+) h(u-) 



2 _ (a(w + ) m(w + ) - b(u+) c(u+)) (h(u_) l(u_) - f(u^) r(u_)) 
a ^ U ' V >- a(u + ) e(u.) e(u + ) h(u.) (CU1J 



3/ \ (- h(u-)) + j{uJ) q(u_)) (e(u+) q(u+) - f(u+) r(u+)) 



, c(«_) i(u_) a(w_)c(w + )e(w + ) (- (^(u_) i(u_)) +j'(u_) m(u_)) 



h{u^)h(u + ) j{u_) a(u + ) e(u-) h(u-) h(u + ) j(u-) 

! (jjOjKj (a(-»+)6(^-)e(^_)r(M + ) -a(M_)&(M + )e(M + )r(M_)) 

a(u+) e(u_) e(u+) h{u ) 2 h(u+) j{u ) 

( , sw x /W2^ . M«_W„^\ W9„^ f &(2u)6( U+ ) , 6("-)(&(«+) 2 -a(«+) <("+)) 



/wo > »,/ w >\ Ur> \ I b(2u)b(u+) . °\ u -) \ °\ u +) ~ a \ u +) H 



a(u+) e(—u-) a(2v) 

V 

| a(u-)b(u + ) (%-)/(«-) -/(M-)r( M -)) | b(u-)f(u-)r(u + ) \ 
a(u+) e(w_) 2 h(u-) e(u_) e(u+) /i(u_) / 

(c(2u)6(2u)) +a(2M)m(2M) /(u_)r(u+)\ &(2tj)c(a_)m(« + ) 
(6(2 u) 6(2 u)) +a(2«)I(2«) + h{u-) h{u+) ) ~ a{2v) h{u_) h{u + ) 

(pin, \h(-i, \ I b(2v)e(-U-)b(u+) \ ( c{2u) c(u-)c(u+) b(2 u) fju.) r(u + ) \ 

a(u+) e(— u_) 



x 
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al(u,v) = - 



f c(u-) m(u+) \ f 6(2 u) b(u+) 6(u_) (b(u+) 2 - a(u+) l(u+)) \ 
\h(u_)h(u + ) J \a(2u)a(u + ) a(u + ) e(uJ) e(u + ) J 



, - (c(2u)b(2u)) + a(2u)m(2u) /(u_)r(u+)\ 
- {b{2u)b{2u)) + a{2u)l{2u) + h{u^) h{u+) ) 



h(n \ ( c ( 2m ) _|_ c(m_)c(m + ) . b(2u)f(u-)r(u + ) \ 
°y U +) \a(2u) + h(u.)h(u+) + a(2u)h(u_)h(u + )J 



a(u+) 

^3 = f(u-)i(u-)ii(u+) e(uJ) f(u+) (- hju-)) + j(u_) gju-j) 

b(2v)f(u + ) } + /(«_)%+) ^ f -(c(2 M )6(2 M )) + a(2 M )m(2 M ) + /( M _)r(«+) ^ 



a(2v)e(M + )y e(u+)h(u-)J V - (6(2 w) 6(2 «)) + a(2 w) /(2 w) h(u_) h(u + ) J 
6(2 «)/(«_)(?(«+) 



a(2 1>) h{uJ) h{u + ) 



r/ \ f - (c(2 «) 6(2 tt))+a(2 u) m(2 «) . /(M_)r(ix + ) \ 
, /(«_)g(u+)\ / V^+/ ^ — (6(2 m) 6(2 u))+a(2 u) 2(2u) h{u-)h(u + )J 

"^• V) = -[ k(u_)k(u + ) ) + ^) (ai6) 

(ai8) 



al(«,u) = 



/ nf—ii \ fin \ ( c ( 2 ") _|_ c("-)c(m+) i 6(2«)/(M_)r(u + ) \ 
ay u-) j yu + ) y a(2u) t ~t~ a {2u)h{u-)h{u+) ) 



a(u+) e{—u ) 

(g(u-)j(u-) -d(u-)i(u-)) fi(u+) c(«+)6(«_) (j(u-)m(u-) - gi(u-) i(u_)) /(«+) 



V 



h(u+)j(u-) a(u+) e(«_) h(u+)j(u-) 

f(u+) (j(u-) g(w_) - iipu-)) (6(w + ) 6(m_) rQu_) /(m + ) - q(m + ) e(^_) c(w_) 

a(w + ) e(w_) e(w+) h{u ) 2 h(u+) j(u-) 
- (c(2 u) 6(2 u)) + a(2 u) m(2 u) /(«_) r(u+) \ 
-(6(2M)6(2u)) + a(2u)Z(2M) + h{u_) h{u+) ) 

' 6(2 m) , 6(m_)6(m + ) \ 
a(2«) e(u_)e(u+) y 

a(u+) e(— u_) I ' e ( u -) e ( u +) h(u-) 

b(u + )b(u^) Z(«_)-/(«_)r(«_)) /(u+)\ 



«(-«-)/(«+) + \ + g( M+ ) (-(/(«_) c(«_)) + /(«-)) 



a(w+) e(-u_) 2 e(tt+) h(u-) J 



(C.19) 
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3 (~ Mu-) i(u-)) + j(u-) m(u_)) (g(u+) h(u+) - c(u+) c(u+)) 

a 9 {u,v) = — — — — — — — 

a[u + ) h{U-) n{u + ) j{U-) 

f(u+) (j'(w_)g(u_) -i(u_)ii(u_)) (6(u+) e(u_) c(u+) r(u_) - a(u+) e(u_) m(u+) r(u_)) 



x 



a(w + ) e(u_) e(u+) h{u ) 2 h(u+) j(u-) 
(a(u+) m(u+) - 6(m+) c(m + )) Q(m_) - /(^_) r(?z_)) 

a(u+) e(u_) e(u+) h{u ) 
f -(c(2u)b(2u)) + a(2u)m(2u) f(u_)r(u + ) \ 

V -(6(2«)6(2«)) + a(2«)Z(2u) h(u-)h(u+)J { ' ' 



101 ' ] ~ e(u + ) h(u.) h(u + ) j(u-) [ ] 



° lK v) " %^7(07K) (C ' 22) 



a\(u,v) = 



d(2 u) 
a(2u) 



d(u—) d{u-\-) 
j(u-)j(u+) 



b(2u)g(u-)g 1 (u+) 
a(2u)j(u-)j(u+) 



c(2u)i(u.)n(u + ) \ bM 
a{2u)j(u-)j(u + ) ) V +' 



X 



6(2 u) b(u+) 
a(2 u) a(u+) 



+ 



a(u+) 

b{u ) (b(u + ) 2 — a(u+) l(u+)) \ 



a(u+) e(u-) e(u+) 





i{u- 


)ii(u+) 


M^(2u) 




M^(2u) 


-)j(u+) 


f 


-)j(u + ) 


M^(2u) 


- + - 





i(u_)ii(u + ) 
j(u-)j(u+) 



+ 



M$(2u) 



M^(2u) 



( c(u-) m(u+) \ 
h{uJ) h{u + ) J 



_ f b{2u)b{u + ) + b{uJ) (b{u + ) 2 -a{u + )l{u + )) \ 
\a{2u)a(u + ) a(u + ) e(u-) e(u+) J 

h(n \ ( c ( 2 ") _|_ c(m_)c(m + ) . b(2u) f{u-)r(u+) \ 
U \ a +) \ a(2u) ^ h(u-)h(u+) ^ a(2u)h(u-)h(u+) } 



M^(2u] 



M$(2u) 



i(2u) h(u-) h(u+) 

a(u + ) 



n(u+) d{u ) \ 
j{u-)j{u + ) ) 



/(«_) r(u+) \ 
h(u-) h{u + ) \ 



(C.23) 
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a\(u,v) 



h(u_)i(u + ) (- {^^f^) +f(u-)H(u + )) h(2v) fi 



h(u+)j(u-) j(u+) 



) fi(u+)g(u-) 
a(2v)j(u-)j(u+) 



| e{u-)h{u_)g{u + )h(u + ) _ H{uJ)h{u + ) _ b(2v) f(u + ) \ 
e{u + )h{uJ) j(u_) j(u + ) \e(u+)h(u_) a(2v)e(u + )) 



x 



i(u-) i\{u + ) 
j(u-)j(u+) 



+ 



+ 



i(u-)i 1 (u+) 


M^(2u) 


- + - 


M$(2u) 


j(u-)j(u^ 


-) 


M${2u) 




M$(2u) 


M^(2u) 




ff(u-)i( 


u-)ii(u+) 



M$(2u) 



h( u -) h(u+)j(u-) 



e(u_) /(«+) (- + g(«_)) 6(2 v) /(«_) g( u+ ) 



e(u+) /i(m_) 2 j(u_) a(2u) %_) %+) 

| /(«_)%+) ^ f -(c(2tt)6(2tt)) + a(2tt)m(2?z) | /(^_)r(^ + ) ^ 



e(u + )h(u_)J V -(b(2u)b(2u)) + a(2u)l(2u) h(u_) h(u+) J 



b(2v)f(u + ) \ 
a(2 v) e{u + ) J 

(C.24) 



4/ \ / /i( M +)#( M -, , 
a£(tt,f) = - — - — — - — - + 



/( u+ ) ( + »m(» i )|mW( 2m) | ^ M<+>(2- i 



j(u-)j(u+) 



+ 



MW(2«) 



M«(2«) 



h{uJ) h{u + ) ) 



e(u+) 
/(«+) 



M$(2u) 



Af<+>(2u) 



Af<+>(2u) 



+ 



/Q_)rQ + ) 
h(u-)h(u+) J 



e(u+) 



(C.25) 



a|(u,u) 



h(u-) h\(u ) i(u+) b{2v) hi(u + ) i(u-) 



x 



M<J(2u) 



a(2u) j{u-)j{u + ) 



M^(2u) 



6(2 u) ^ + i(w_)j(w+) ^ 



a(2u) /i(w+) / h(u + ) j(u-) J 

(C.26) 



a 7 (w, f ) 



j( u -)j( u +) J 



i(u+) 



i(u-) h(u + ) 
j(u-)j(u+) 



h(u+) 



M<+ ) (2«) 



M<+ ) (2n) 



(C.27) 
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<4(u,v) 



n{uS) {a(u + ) bi(u + ) — d(u + ) c(u + )) 
a(u+)j(u-)j(u + ) 



+ 



i(u-) ii(u+) 


M$(2u) 




M$(2u) 



+ 



x 



M^{2u) 

(o(tt+) m(u+) - b(u+) c(u + )) (fe(tt-) IjuJ) - f(uJ) r{u_)) 

a(u + ) e(u-) e(u + ) h{u_) 
f\(u_) g(u+) (- (a(u+) ejuJ) m(u + ) r{uJ)) + b(u+) e{uJ) c(u+) r(u_)) 

a(u+) e(u-) e(u+) h{u )j{uJ) j(u+) 



- h 1 (u-)i(u+) 



gi(tt-) (a(u+) h(u + )-c(u + ) c(u + )) _|_ f(u+) ii(u-) (b(u+) e(u-)c(u+) r(u-)-a(u+) e(u-) m(u+) r(u-)) 



X- 



a(u+)h(u+)j(u-) 



a(u+) e(u_) e(u+) h(u~) h(u+) j(u-) 



j{u_)j{u + ) 



i(u-)i 1 (u + ) 
j(u-)j(u+) 



^-( gi (u^)i(u-)) + j{uJ)m{uJ)) (a(u+) h(u+) - c(u + ) c(u + )) 
a(u+) h(u-) h(u+) j(u-) 
f(u+) (j(u-) q(u-) — i(u ) ii(u-)) (b(u+) e(u-) c(u+) r{u ) — a(u+) e(u-) m(u+) r(u-)) 



a(u+) e(u-) e(u+) h(u-) h(u+)j(u 
{a(u + ) m(u+) - 6(u+) c(u+)) (/i(u_) Z(u_) - /(u_) r(it_)) 



M<+>(2u) 



Af<+>(2u) 



/(«_)r(«+) \ \ 
h(u-) h(u+) j 



a(u + ) e(u-) e(u+) h(u-) 



(C.28) 



a?o( w > u ) 



j( M -) 2 i( M +) 



+ 



j(u-)j(u+) 



h{u-) g{u+) gi(u + ) I i(u_)ix(« + ) 



<(2«) 



x 



(e(tf+) - ffi(^+)) (j(u_) q(u_) - i{u_) h{uJ)) 

e(u+) h(u-) h(u+) j(u-) 

a(— U-) h(u+) 



b\(u,v) = 



a(u + ) h(—u_) 



(C.29) 



(C.30) 
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bl(u, v) = (e(u+) q(u+) - r(u+) f(u+)) 

{u-)-g(u-)q(u-)) \ j_ -(g(u-) 

e(u~) i{u-) 



X 



(-(91 ( u - ) ))+.?'("- ) )) ) Mu- )-g(tt- ) g(tt_ )) \ , m(u_))+t(u_) g(tt_) 

e(u_)i(u_) (-(i(u_)ii(«_))+j(u_)g(u_)) y 



e(w + )e(M+) 

(C.31) 

bl(u,v) = (/ii(u_)Zi(u_) - /i(«_)ri(u_)) 

/ _ m(u + ) (a(u+) m(tt + )-6(u + ) c(n + )) _ c(ix + ) (-(b(« + ) m(u + ))+l(u + ) c(u + )) \ 

\h(u+) h(u + )(~(b(u+)b(u + ))+a(u + )l(u+)) h(u + )(-(b(u+)b(u + ))+a(u+) l(u+)) J fnQ0 \ 

X h{uJ)hx{uJ) (C - 32) 

e (u ) ( ei ( u+S) — g( M +)(-(gi(^+)g("+))+/i(^+)^("+)) _ (e("+) /1 ("+)-/("+) gi("+)) n(«+) \ 

IV ^ j(u+)(-(e(u + )g(rt + ))+/(tt + )r(w + )) j(it+) (e(tt+) q(u+)-f(u+) r(tt+)) 7 
«) = ^ ^- 33) 



1 / \ A 



(C.34) 



c\Au,v) = - ( b f u \ b M ) _ (c.35) 
12V ' ; \a(2v) a(u+)J a(u + )e(u_) v ; 

U(n \ ( _ ( b(2v)b(u + ) \ , b(u-)e(u+) \ 
= °y ZU > \ \a(2v)a{u + )) + g(n + )e(n_) y ) b(u+) e{uJ) _ b( U ^) b(2 v) e(u + ) 

l[U,V) a(2u) a{u + )e{u_) a(2 v) a(u+) e(u_) 1 J 

"(S) (c ' 38) 

)/("-)&("+) _ /V,. \ If,. \ _ /(n+) (-(c(2 u) 6(2 M ))+a(2 n) m(2 u)) 

„2 /, , , x _ a(«+) il M -JH^+J -(h(2 U )6(2tt))+a(2u)t(2u) ( „ . 

° 22 ^ Uj " e(u + )h(u.) ( j 

cl lM = (- f^)fM-)) + ( b ^\ b ^\ + /(u _ )6(t|+) 



a(2«) y ya(2«)a(« + ) a(u + )e(u-) 

b{2u)f{u^)l{u + ) 



a{2u) 



- c(u_) /(«+) /(e(« + )M«-)) (C40) 
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1.3 



_ f b(2v)f(u + ) \ f(u.)h(u + ) 
1 ' ' " \a{2 v) e(u+) / e(«+) 



(C.41) 



c? 2 (u,u) 



6(2 
a(2 



e(u + ) h(u-) 



e(u+) h{u ) 



/(u_) /i(u+) 




6(2 u) 


+ 

e(u+) h{u ) 


Mg } (2 M ) 





/(u + )M«-)|M(+'(2«) 



itf<+>(2u) 



a(2 i>) e(u+) 



(C.42) 



c(2ti)/(^)% + ) 
a(2«) 



°\, Z "M la(2») a ( U+ )J t a(«+)e(«-) y ) ; e(li_) 

a{u + ) e{u) 



a(2u) 



+ 



a(2 f ) a(u+) e(u_) / 
/("+) 

o(2«) 



^2 -y) ( c ( 2 ")/("+)M"-) _|_ b(2 «)/(«-)/("+) 



a(2u) 

,,,2,,, +c(u_)/(«+) 



a(2v) 



/(e(u + ) /*(«_)) 



(C.43) 



C2 4 (m,w) = - 



(C.44) 



4»( u > u ) = - 



*(«+) 


M<J>(2u) 


?(w_)g(w+) 


h(u+) 


<'(2u) 


h(u + )j(u_) 



+ 



(C.45) 



5(m + )^(m_)c(m + ) 
a(u+) h(u + )j(u_) 



i(u + )j(u-)M+ ) (2u)\ i(u-)\M^(2u)\q{u + ) 

g(u_) m{u + ) H i .. j 1 + 



M<+>(2u) 



M«(2n; 



h(u + )j(uJ) 



\ I b(u + )f(u-)b(u+) ,/ x,/ x f{u+)h{u-){-{c{2u)b{2u))+a{2u)m{2u)) \ / x 
H M -J I ^) A«-JH M +J -(b(2u)6(2 U ))+a(2«)J(2i i ) ) r \ U +) 

e(u+) h{uJ) h{u + )j{uJ) 



(C.46) 



47 



c 21 (u,v) 



b % u \ b{ ? + \ + b( "-j e( " + ^(u_) c(u + ) 

a(2u) a(u + ) a(u + )e(u_) / ^ v / V tv 



h(u+)j(u-) 

d(2u)i(u^) j(w-) _|_ 6(2«)j(u_)m(ii+) j_ c(2 u) ) 



a(2u) 



a(2 u) 



+ 



a(2 u) 



r(-u + ) 

6(2^)/(^-)f(^ + ) 
a(2 u) 



h(u + )j(u-) 
c(2u)f(u + )h(u.) 
a(2u) 



a(2 u) a(u+) a(w + ) e(w_) / 



-c(u_) /(«+)) / (e(u + )/l(u_)/l(u+).7'(u_)) 



11 



_ / b(2v)i(u+) \ i(u-)j(u+) 
1 ' ] ~ \a(2v)h(u + )J + h(u + )j(u.) 



(C.47) 
(C.48) 



[u,v) 



i(u-)j(u+) 


M$(2u) 


h(u+)j(u-) 


<(2«) 



+ 



i(u+)q(u-) 



6(2 v) 



h(u+)|M(+ ) (2u)| 



<?("+) 
/i(tt + )j(u_) 



h(u + )j(u-) 

( ( b(2v)f(u+) \ f(u-)h(u+) \ ., x / x 
^ ^a(2^)e(« + ) J + e(u+)h(ti-) J H u -J r l M +J 

h(u + ) j(u-) 
g{u + )m{u_) 



a(2v) 



(C.49) 



C 12 [U,V) 



+ 



i(w_)j'(u+) 


M$(2u) 


h(u + )j(uJ) 


M$(2u) 



( (b(2v)b(u+)\ 
\ \a(2v)a{u+)) 


1 £>(«-) e 0+)^ 
I" a( U+ )e( U _)J 


g(u_)c(u+) i(u+) 


M$(2u) 


q(u-) 


h(u+)j(u-) 


h(u + )j{u_) 


M$(2u) 



b ( 2v ) ( 9(n-)m(u + ) i(u+)\M { + i \2u)\ i(u-) \M%>(2u)\ g{u + ) 
\ > yh(u+)j(u-) |m : <+ ) (2m) |M^+ ) (2n) 



r (+), 



a(2u) 



- Mm 



- ( b ?A K ? + \ ) + txN) /(«-) Ku+) 



e(u+) h{u ) 

/(«_) %+) (- (c(2 u) 6(2 u)) + a(2 u) m(2 u)) 
e(u+) h{uJ) (- (6(2 u) 6(2 «)) + a(2 w) Z(2 «)) 

6(2.) (/(«_) f( M+ ) + /( ^%jjff4tgl^ g^) 
a(2 f ) e(u+) h{u ) 



+ 



/(«+)*(«-) 
e(u + ) h(u-) 



r(u+) ] / (h(u+)j(u-) ) 
(C.50) 
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3 / \ d{u + )g(u_) d{2u)i{u_)j(u + ) 6(2 u) g(u+) m(u_) 

C ( W ) — — i- — i— 

h(u + )j(u-) a(2u) h{u + ) j{u^) a(2u)h(u + ) j{u ) 

( b(2v)b(u+) \ , b(u_)e(u + ) \ \ 

U(2"M"+)J + °("+M"-)J , 6(u+)e(u-) _ fc(«_)fc(2t))e( M+ ) I / \ / s 
a(2u) a(«+)e(u_) a(2 u) a(« + ) e(u_) J J\ —J \ +/ 



...... h/'OiA ( ^«)'("+) I | t(2«)g(n-)m(u + ) , c(2u)i(u-)q(u+) \ 

C(2u)l(u + ) q{U-) y a (2u)h(u + ) ~T h(u+)j(u-) ^ o(2«)/i(ti+)i(«_) ^ o(2u)A(«+)j(«_) / 



a{2u) h{u + ) j{u) a{2v) 

I h(0 \ ( — ( b ( 2 ^) b ( u +) \ I b(M~)eQ + ) 

C(2tt)/(«_)|j(tt + ) / °l Z-U J ^ \a(2v)a(u+)J + o(u+)e(u_) 



a(2tt) I a(2tt) 



6(w+)e(w_) 6(w_)6(2t;)e(w + )^ 6(2 tz) /(«+) Z(m_) , , , , 

}{u-)b{u + ) H — hc(« + )/(u_) 



a(u+) e(w_) a(2 v) a(u + ) e(u-) J a(2u) 
b(2v) ( «*»)f<£)««-) + ^j^ptl + c(u _) 



a(2u) 




/ ( e(«+) /i(w_) /i(«+) j( M -)) 

(C.51) 



*->~(S) 



(C.52) 



(-(/( M _) C (n_)) + %-)/(«-)) \ 6(2 «)/(«+) 
V e(«+) (/i(u_)Z(u_) -/(«_) r(u_)) / a(2u)e(u+) 



(C.53) 



dl(u,v) 



a(u+) m(u+) — b(u + ) c(u+) 
-(b(u + )b(u + )) + a(u + )l(u + ) 



(C.54) 



M$(2u) 


(a( 


u+) m(u + ) 


- 6(u+)c(«+)) 




M$(2u) 





d|(u,u) 



(- (^i(«-) *(«-)) + j( M -) m ( M -)) 5(«+) , &(«+) r(u-) f(u+) 



- (i(u_) ii(u_)) + j(u_) g(u_) 



+ 



/i(w_ 



/ a(« + ) r (rt+ ) (-(gi (rt_ ) ))+j(u- ) to(m_ )) _ b(u+) e(u+) r{u-) \ ft \ 
\ -{i(u-)h{u-))+j(u-)q(u-) h(u-) ) J\ a +) 



e(u + ) 



I (a(u + )l(u + ) -b(u + )b(u + )) 



(C.55) 
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W">-\ a (2u) 
a(u+) ii(u+) (- (d(u-) i(u-)) + gju-) j(u-)) 
- (i(u_) ii(u_)) + ?(«-) 

a(2u) (-(i(u_)ii(u_)) + j'(u_) g(u_)) 
h{u + ) (— {l{u) c(u_)) + r(u_) f(u-)) 



+ »n(u_) - 0i (u_) 



x 



-(%-)*(«-)) + /(«-M«-) 

' \^ e(u + )(/i(u-)Z(u-)-/(u_)r(w-)) r a(2 U )e(« + ); 

- (i(w_) + j'(u_) g(u_) 

/(-(&(«+)&(«+)) +a(« + ) /(«+)) (C.56) 



= - I ,K)^)-/WrK) J (a57) 



d|(u,u) 



(-(e(u + )/ 1 (u + )) + /(u + )fc(u + )) 


M«(2«) 









+ e(u+) /i(u+) 



/ li(tt+) m(-» + ) (a(-»+) m(tt + ) - b(u + ) c(u+)) 
U(«+) (" (&(«+) &(«+)) + *(«+)) 

c(« + ) (-(&(«+) m^ + Z^cOOA (fci(u _ )Zi(li _ ) _ /i(li _ )ri(u _ )) 



(- (6(«+) &(«+)) + a(«+) Z(u+)) ; 
/ (/n(«_) (/!(«_) n(«_)-/n(«_) Zx(«_))) ) / (e(« + ) (?(«+) -/(« + )r(« + )) (C.58) 
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e(«+) fi(u-)k(u+) 


M$(2u) 


(e(u+)/i(«+) - 

i 


/(u+)^i(«+)) 


M«(2 M ) 


hi(u-) 


M$>(2u) 




M 2 ( J(2 M ) 





/(u+) (ii(u_)m(u_) -£i(u_)g(u_)) g(u+) , ra(u_) e(u+) /i(u+) 



+ 



/ rt(M-)e(M + )gi(u + ) \ . f(u+)r(u+) (n («_) m(«_ )-gi (u- ) \ f / \ 

^ /n(«_) J i(«_)u(u_)-j(u_) g (tt_) J JK U +) 



e(u+) 



M^\2u) 


(a(u + ) m(u + ) 


-b(u+)c(u+)) 




M$(2u) 





e(u+) f 1 (u-)m(u+) 



a ( u +) (- (9i( u -) *(«-)) + j( u -) m(u-)) q{u + ) b(u + ) r(«_) /(u+) 



(u + )r(u + ) (-(gi(u-) t(M-))+j(«-) m(tt-)) 6(m+) ep + ) rQ_) 



7i(u_) 



h(u. 

/(«+)' 



e(u+) 

/ (/ii(u_) (a(u+) Z(u+) - 6(u+) 6(u+) ) ) + c(u+) e(u+) /i(u_) 



x 



M 2 ( J(2 M ) 


(&(«+) m(u+) 


-Z(«+)c(u+)) 




M<J } (2u) 





6 ( u +) (- *(«-)) + m(u_)) g(u+) Z(«+) r(«_) /(u+) 



-(i(u_)ii(u_)) h{u_ 

f r{u+)b{u+) (-( gi (u-)i(u-))+j(u-)m(u-)) e(u+) l(u+) rju-) \ ft \ 
\ -(i(u-)h(u-))+j(u-)q(u-) h(u-) ) J\ a +) 



e(u+) 

I (/ii(u_) (a(u+) Z(u+) - &(u+) &(u+)) ) ) / ( e(u+) g(u+) - /(«+) r(u+) ) 



(C.59) 
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di(u,v) 



d(2u) (- (e(u+) /i(u+)) + f(u + )g 1 (u+)) 



a(2u) 

f(u+) h(u + ) jjg{uJ) ijju-) - d(uJ) g(u-)) 
i(u_)ii(u_) - 

g(- M _) e( M+ ) A( M _) %+) ( gggj^ + c(iQ c( M+ ) + K2u) %i riu+) ) 
a(u+) h(-u-) (- (/ii(u-) ii(u-)) + 

^/g^N f /("+) m(ix_)-gi(u_) qp + ) _ ep + ) ) ri(«_ )-n(u_) h (u_ )) /i \ 

v ' y i{u-) h(u-)—j(u-)q(u-) -(hi(u-)h(u-))+fi(u-)ri(u-) J 



a(2u) 

_ ( f(u+) r(u+) (ii(u-) m{u-) - gi(u-) g(u-)) 

V i(u_)*i(«-) -j(u-)q(u-) 
e(u+) gi(u+) (fi(u-) nju-) - n(u-) h(u-)) 

- (/n(u_) Zi(u_)) + n(u_) 

(-(/(tQ c ( M _)) + ft( M _)/( M _)) 6(2tQ/( M+ ) \ 

V e(u + ) (h(u_) l(u_) -/(«_) r(«_)) a(2 u ) e( u+ ) J 
/(e(« + )9(« + )-/(« + )r(« + )) (C.60) 
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4 fb{2v) n{u + ) d(u-)\ a{uJ) n{uJ) d(u + ) e{u + ) 

2 ' \a(2v) j(u_) j(u + )J a(u + )e(u-)j(u-)j(u + ) 

( d{2u) _(__) d{u + ) b(2u)g(u_)g 1 (u+) c(2u)i(u-)i 1 (u + ) \ 
\a(2u) j(u-)j(u+) a(2u)j(u-)j(u+) a(2 u) j (u_) j (u+) J 
^ e{u + )b{-u_) | b{2v)b{u + )\ 



_(_+) e(— __) ' a(2v)a(u+)J 
/_(--) _?(-+) (a(-+) 6(m_) e(__) r{u + ) - ajuJ) b(u + ) e(u + ) r(__)) 
a(u + ) e(__) e(_+) h{uJ) j(u-) j(u+) 

M$(2u) 

+ - 





i(u-)i 1 (u+) 


M^(2u) 


-)j(u + ) 


f 

j(u-)j(u+) 


M$(2u) 



( < \U \ fb(2u) , &(__)&(_+ A b(2v) ( b ( 2 ^ b ^ b(u-) (b(u + ) 2 -a(u+)l(u + )) 

(_ + )6(-__) ( jaliy + 4 U _) e(u +) J ^^UP")^) »(«+)_(«_)_(«+) 



X 



V 



a(_+) e(— __) 



+ 



a(2u) 



o(h-)6(m + ) (%,)l(^)-/(M_)r(M.)) 6(n_)/(n_)r(n + ) > \ 
+ a( M+ )e( M _) 2 %_) + _(__)_(_+)*(«-),/ 



/(-+) u(«-) (a(«+) 6(u_) e(u_) r(u+) — a(u_) 6(u+) e(u+) r(«_)) 
a(w+) e(u_) e(u+) h(u~) 



i(u-) ii(u+) 
j(u-)j(u+) 



+ 



M^(2u) 



M^(2u) 



h(u + )j(u-) 2 j(u+) 

c{u ) gi(u + ) i(u-) b{2 v) c(__) m(u+) 
h(u-) h(u+) j(u-) a(2v) h(u-) h(u+) 



o(tt_) c(u+) e(u + ) {j{u_) m(__) - gi{u_) i{uJ)) 

a{u + ) e(__) h{u_) h(u + ) j{u ) 
fX u +) U( u ~) q( u -) - ( a ( u +) K u -) e ( u ~) r ( u +) - a ( u ~) K u +) e ( u +) r ( u -)) 



a(u+) e(u-) e(u+) h(u-) h(u+)j(u-) 



( 



i v ) h (-v ) bi2u ^ i b{u-)b{ U+ ) \ b(2v) y 9 u t + 



(«-)e(w+) 7 



V 



a(_+) e(— 



+ 



b(2 -)6(_ + ) i b ( u ~) (K u +) 2 - a ( u +) l ( u +)) \ 
| ______ _____ i 



a(w-|-) e(u_) e(«+) 



a(2v) 



a(-_)6(-+) f(-_) - /(__)r(__)) &(__) /(«_) r(u+) 



a(u+) e(u-Y h(u-) ' e(u-) e(u+) h(u-) 

(c{2u)b{2u)) + a(2u)m(2u) /(u_)r(_+)\ 
X ' -{b(2u)b(2u)) + a(2u)l(2u) + h{uJ) h{u + ) ) 

An, \h(-i, \ I b(2v)e{-u-)b{u+) \ ( c(2u) , c(u-) c(u+) , b(2 u) /(__) r(u+) \ 
■\ U +) U \ a ~) a(2«) y ^a(2u) ^ /i(«_)h(_+) ^ o(2_) h(__) h(u+) ) 

a(u + ) e(— __) 



(C.61) 
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at(u,v) = - 



) , d(u-)d(u+) _|_ b(2 n)g(u-)gi(u + ) , c(2u)i(«^)ii(u + ) \ 
) a(2u)j(u_) j(u+) o(2«)i(ii-)i(ii + ) y 



a(u+) e(— w_) 

ff^+j (6(m + ) 6(m_) r(u_) /(«+) - a(u+) e(^_) c(m_) g(u+)) 
a(«+) e(u_) e(«+) h{u_) j{uJ) j(u+) 



n( 



a(u+)e(u-)j(u-)j(u+) 



x 



,r_ 7/ wr 7/ n f 6(2 «) , b(«-)fc(»+) A 

K M -J / L(2«) + e(u_)e(u+) J 



+ 



+ 



i(u_)ii(u+) 


M$(2u) 




M$(2u) 


J'(w-) 3i u +) 




- + - 


M^(2u) 



a[u + ) e{—u^j 



g{u + ) (-(f( u _)c{u_)) + h{u_)f{u_)) 
e(u-) e(u+) h{u ) 



b(u + )b(u_) (h(u_)l(u_)-f(u_)r(u_)) f(u + ) \ 
a{u + ) e(uJ) e(u+) h{u ) ) 

d(u-)fi(u+) c(u+) gi («_)&(«_) /(«+) _| f(u+) i(u- ) (b(u+) b(u-) r(«_) f(u + )-a(u+) e(«_) c(u-) <?(«+)) 



X 



h(u+)j(u-) a(u + )e(u-)h(u+)j(u- 



+ 



a(u+) e(«_) e(«+) h(u~) h(u+) j(u~) 



i(u-)ii(u+) 
j(u-)j(u+) 



+ 



M^(2u) 



M^(2u) 



- (d{u-)i{u-)) + g{u-) j{u-)) f 1 {u+) 
h{uJ) h(u+)j(u-) 



a —u 



\ fir,, \ ( c ( 2 ") _|_ c{u-)c{u+) 6(2 u)f(u-)r(u+) \ 
-) J \ a +) y a (2u) ^ h(u-)h(u+) ^ a(2u)h(u-)h{u+) ) 

a(u+) e(— u ) 
c(u+)b(u-) (- (gxiu-) i(u-)) +j(u-)m(u-)) f(u+) 



+ 



a(u+) e(u-) h{u ) h(u+) j(u-) 
f(u+) (j(u-) g(w_) - i{u_) ii{u-)) (b(u+) b(u_) r(u_) f(u + ) - a(u + ) e{u_) c(uJ) g(u + )) 

a(u+) e(u-) e(u+) h{uJ) 2 h(u+) j{u ) 
- (c(2 u) b{2 u)) + a(2 u) m(2 u) /(«_) r(u+) \ 
-(&(2u)&(2u)) + a(2u)Z(2u) + /i(u_) /i(u+) J 



x 



6(2 u) b(u_)6(u + ) > \ 



a(2u) + e(u-)e(u+) M g(« +) _) / ' (u _) - / "(u _) c(u _)) 



a(w+) e(— 



+ 



e(w_) e(u+) h{uJ) 



b(u + )b(u-) (h(u-)l(u-)-f(u-)r(u-)) f(u + ) \ 



a(u + ) e(/u_) e{u + ) h(u. 



(C.62) 
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Appendix D: Relations for arbitrary S 



In this Appendix we present certain expressions concerning the unwanted terms of the one- 
particle problem as well as the construction of the two-particle vector for arbitrary S. 

The commutation rules used in the solution of one-particle eigenvalue problem come from 
the entries [1,2], [2,3], [2S,2S+1], [2,2+2S+l], [3,3+2S+l], [2S+1,2(2S+1)] of the 
boundary Yang-Baxter equation (jB.lj) . To cancel the unwanted terms we need to know how 
to compute the ratio %y ^ A:A 



flF(A,Ai) 



which is not expected to have a dependence on the i-th index. 



This means that this ratio can be calculated collecting the simplest unwanted contributions 
which turns out to be those coming from the commutation rules between the fields A 2 s{ty, 
^2S+i(A) and -Bi 2 (Ai). Considering the help of mathematical induction we find that the func- 
tion q^g(X, ^1) 1S 

. -"-1,25 



qfg (A, Ai 



-uj 



(+) 

2,9 



(A)- 



? 25 + l ,2 (A + A 
-^l,2s(^ + ^1 



l > (+) 



( 



p2S+l,2 
•"-1,25 



(A + A x ) 



M (+) 

IV1 2S,2S+1 



(2A) 



d2S+1,2 
-"■1,25 



(A - Ai) 



^2,2S+l(^ + ^l) ^2,2S+l(^ + ^l) 



R 1,2S(^ + ^) R h2S+l(^ + ^1 



25+1,1, 



V 



/ 

(d.i; 
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while (A, Ai) is given by 
gg(A,A 1 ) = 4+ ) (A) 



,25+1,2, 



p25+l,l, 
-"■1,25+1 ^ 



^l) ^1,25+1 + ^0 


^1,25 1>2 (^ + ^0 


M#(2Ai) 


- Ai)i2i 2 ' 5 (A + Ai) 


R l,2s(^ + ^l) 


M 1 ( + ) (2A 1 ) 



(+) 

25+1 



(A) 



M 



(+) 

25,25- 



-i(2A) 



,25+1,2/ 



m£L(2A) 



,25+1,1, 



^i)-^i,25+i(^ + -vo 


^l,25 1,2 (^ + ^0 


M 1 ( + ) (2A 1 ) 


- Ai)i?i 2 ' g (A + Ai) 


R 1,2S(^ + ^0 


M 1 ( + ) (2Ai) 



M 1 (+ 2 ) (2A 1 ) 


^1,25 1,2 (^ ~~ ^l) 


R i,2s(^ + ^i) R i%; 1,2 (^ + ^l) 

^2,25+l(^ + ^l) -^2^5+1(^ + ^1) 


M 1 ( + ) (2A 1 ) 


^l,25+l(^ ~~ ^l) R \,2s{^ + ^l)-^l,25+l(^ + ^l) 



~ ^l)R%2S i,Z (^ + ^l)-^2^2S+l(^ ~ ^l) _ ^iS^i^ ~~ ^l)-^25+l(^ + ^l) ^%2S^ ~ ^1 



,25+1,2, 



,25+1,2 



25+1,2, 



,25,1 



,25+1,2, 



+ 



^l,25+l(^ ~ ^l)-^l5s(^ + - ^1 



,25,1 



25+1,1, 



(D.2) 



Taking into account the explicit expressions for the Boltzmann weights and for the functions 

(2), 
2S \ 

g^(A,Ai) 

We close this Appendix discussing the construction of the two-particle state. The appro- 



ojf(X) one can verify that the ratio 7rj7~^ satisfy Eq. (|107jl . 



priate commutation relation is derived by combining the entries [1,3] and [1,3+2S] of Eq. ljB.lj) . 
After some algebra we find that commutation relation between the creation operators Bi 2 (u) 
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and Bi 2 (v) is 



Rl,2{ U +) ~i R 1,2( U +) 

25+1 pi, 2 / \ 



B l2 {u 



1,2 ( 



-12O) 



= Z s (u,v) B 12 (v)B 12 (u) + — — 



+ 



Jig(«_) J2jJ(u_) 

#2,2 (W-) #2,3 ( W - 




,1,2/ \ 25+1 D l,2 / \ 

#2l( M +) ~i #2l( M +) 



^3l(%) „ / v. 

j R 21 (m+) 



25+1 t,1,3 / \ 
j =4 ^2,ll W +J 



21 1 U 



(D.3) 



where function Zg(u,v) has been defined in Eq. (jll3|) . 

The above relation allows us to define the following vector 

Si2(«)Si 2 (v) + -2^7 — -5 13 (m)A 2 (v) + 2^ — : 

#1, 2 (W+) ~i Rl,2{U+, 

Rl' 2 2 (u_) (r*>*(u+) 2 ^ R{%-\u + ) 

1 -Bi 3 (u)Ai(v) + —pjTJT- — ~Bij(u)Cj-. 2 i(v) ) (D.4) 



[u,v] 



-Bij{u)Cj- 12 {v) 



Ri, 3 (u-) \Rt^ 



RMu + ) 



which is symmetric under the exchange of the variables u and v, thanks to certain identities 
between the Boltzmann weights. More precisely, we have 



(u,v) = Z s (u,v)(f)(v,u) 



(D.5) 



The two-particle state is now obtained by acting the vector (jD.4j) on the pseudovacuum 
|0s) leading us to 



|^2(Ai, A 2 )> 



j?i2(A,)/.j, : .(,Ao) f + ^ B v ,{ A, IU ,Ao) 



#i^(Ai + A2 



R^(X t - A 2 ) ^(A x + A 



3,1, 



#1^3 (Ai — A 2 ) -Ri^ 2 (Ai + A 2 
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-S 13 (A 1 )A 1 (A 2 ) |0 5 > 



(D.6) 



Finally, taking into account Eq.(|90|) we then recover the expression exhibited in 

section 14.31 

References 

[I] P.P. Kulish, N.Y. Reshetikhin and E.K. Sklyanin, Lett. Math. Phys. 5 (1981) 393 

[2] H.M. Babujian, Nucl.Phys. B, 215 (1983) 317; LA. Takhtajan, Phys. Lett. A 87 (1982) 
479 

[3] K. Sogo, Y. Akutsu and T. Abe, Prog. Theor. Phys. 70 (1983) 730 

[4] R.J. Baxter, "Exactly Solved Models in Statistical Mechanics", Academic Press, New York, 
1982. 

[5] V.E. Korepin, G. Izergin and N.M. Bogoliubov, "Quantum Inverse Scattering Method and 
Correlation Functions", Cambridge University Press, 1993 

[6] I. Cherednik, Theor. Math. Phys. 61 (1984) 977 

[7] E.K. Sklyanin, J. Phys. A: Math.Gen. 21 (1988) 2375 

[8] L. Mezincescu and R.I. Nepomechie, J. Phys. A: Math.Gen. 24 (1991) 211 

[9] L. Mezincescu, R.I. Nepomechie and V. Rittenberg, Phys. Lett. A 147 (1990) 70; CM. Yung 
and M.T. Bachelor, Nucl.Phys. B, 435 (1995) 430 

[10] A. Doikou, Nucl.Phys. B, 634 (2002) 591 

[II] R.I. Nepomechie, J.Phys.A.Math.Gen. 37 (2004) 433; J.Stat.Phys. Ill (2003) 1363 
[12] J. Cao, H.Q. Lin, K.J. Shi and Y. Wang, Nucl.Phys. B, 663 (2003) 487 

[13] H. Fan, BY. Hou, K.J. Shi and Z.X. Yang, Nucl.Phys. B, 478 (1996) 723 



58 



[14] R.I. Nepomechie and F. Ravanini, J.Phys.A: Math.Gen. 36 (2003) 11391 
[15] J. de Gier and P. Pyatov, J.Stat. 03 (2004) P002 

[16] H.J. de Vega and A. Gonzalez-Ruiz, Mod.Phys.Lett.A, 9 (1994) 2207; J.Phys.A: 
Math.Gen. 26 (1993) L519 

[17] T. Inami, S. Odabe and Y.Z. Zhang, Nucl. Phys.B, 470 (1996) 419; A. Lima-Santos, 
Nucl.Phys.B 644 (2002) 568 

[18] L. Mezincescu and R.I. Nepomechie, J.Phys.A: Math.Gen. 25 (1992) 2533 

[19] G.A.P. Ribeiro and M.J. Martins, pimST/0406021 , Nucl.Phys.B (2004) in press 



[20] R.B. Stinchcombe and G.M. Schutz, Phys. Rev. Lett 75 (1995) 140 

[21] M.T. Batchelor, in Proceedings of the 22 n d International Colloquium on Group Theoretical 
Methods in Physics, Eds. S.P. Corney et ai, International Press, Boston (1999) 261. 

[22] N. Kitanine, J.M. Maillet and V. Terras, Nucl. Phys.B, 554 (1999) 647 

[23] Y.S. Wang, Nucl. Phys.B, 622 (2002) 633 

[24] H. Fan, Nucl.Phys.B 488 (1997) 409 

[25] V.O. Tarasov, Teor. Math. Phys. 76 (1988) 793 

[26] M.J. Martins, Nucl.Phys.B 450 (1995) 768; M.J. Martins and P.B. Ramos, Nucl.Phys.B 
500 (1997) 579 

[27] G.L. Li, K.J. Shi and R.H. Shi, Nucl.Phys.B 670 (2003) 401; G.L. Li, K.J. Shi and R.H. 
Yue, Nucl.Phys.B 687 (2004) 220 

[28] X.W. Guan, J.Phys.A: Math.Gen. 33 (2000) 5391; A. Foerster, X.W. Guan, J. Links, I. 
Roditi and H.Q. Zhou, Nucl.Phys.B 596 (2001) 525 



[29] V. Kurak and A. Lima-Santos, Nucl.Phys.B 699 (2004) 595; \nlin.SI/04 07006\ 

59 



[30] A.B. Zamolodchikov and ALB. Zamolodchikov, Ann.Phys. 120 (1979) 253 

[31] S. Ghoshal, Phys.Lett.B 334 (1994) 363; E. Corrigan and Z. Sheng, Int.J.Mod.Phys.A 12 
(1997) 2825; M. Moriconi, Nucl.Phys.B 619 (2001) 396 

[32] P. Fendley and H. Saleur, Nucl.Phys.B 428 (1994) 681; C. Ahn and R.I. Nepomechie, 
Nucl.Phys.B 586 (2000) 611 



[33] W. Galleas and M.J. Martins, nlin.SI/0407027 



60 



Manifold 


e = e+/e_ 




I 


+ 


Upper 


II 




Triangular 


I 




Lower 


II 


+ 


Triangular 



Table 1: The triangular property dependence of K s (A) on the ratio 
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